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BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

This invention relates to ultrasonic acoustic imaging, primarily for medical purposes. 

2. Brief Description of the Background Art 

Ultrasonic acoustic imaging finds many uses, particularly in the field of non-invasive 
medical testing. Direct detection of the emitted acoustic frequency permits, for example, 
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prenatal fetal imaging. Detection of Doppler shifted acoustic frequencies permits observation of 
flow of a particle-containing fluid, such as blood flow. Acoustic imaging equipment utilizes a 
probe that is applied to the skin of the patient overlying the part of the body being investigated. 
At the end of the probe is one or more transducers, usually piezoelectric, that are excited by 
bursts of electrical energy at the ultrasonic frequency and, for example, square wave modulated 
to transmit square wave modulated bursts of ultrasonic energy into the body region being 
investigated. The subsurface structures reflect some of that energy, either at the transmission 
frequency or Doppler shifted, back to the probe, where it is detected by piezoelectric receiver 
elements in the probe. 

One application of this technology to the three dimensional mapping and tracking of 
blood flow is disclosed in parent US Application No. 09/926,666, which is scheduled to issue on 
January 27, 2004 as US Patent No. 6,682,483. The pertinent text of that application is included 
herein below. 



SUMMARY OF THE INVENTION 

When studying blood flow in a subject, the magnitude and phase of the Doppler shift 
detected in the reflected acoustic energy indicates amount of blood flow and the velocity and 
direction of blood flow within a blood vessel. The signals received by each of the detector 
elements in an acoustic probe is analyzed so as to focus the receive beam to examine the 
reflection from each spot as a series of acoustic energy bursts is directed to each spot. The 
received signals are passed through a Wall filter to filter out the unshifted acoustic energy 
reflected by stationary tissue. The magnitude of the autocorrelation function of the reflection 
from each pulse is used to indicate the amount of flow and the phase of the autocorrelation 
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function comparing the signals from one pulse to the next is used to determine the velocity and 
direction of flow. Analysis of the received signals by comparing the signals from one pulse with 
the signals from the next pulse in the burst is termed "analyzing the signals at a lag of one". The 
inventive method disclosed herein uses the magnitude of the autocorrelation function at a lag of 
one as an indication of flow at that point. That information is stored and/or displayed as one 
point in the total analysis of the subject. 

DETAILED DESCRIPTION OF THE INVENTION 

The output of the Wall Filter is used for Power Doppler. If multiple pulses or multiple 
samples are available, the power at a particular location is commonly obtained by averaging the 
squared magnitudes of the samples. If the I th complex sample (amplitude and phase) is denoted 
by 2/ , i = 1, 2, . . ., N 9 the power is 



The phase of the autocorrelation function, a(f) or a T9 evaluated at a lag of one (r= 1 PRI), is 
commonly used for Color Doppler. In this invention, we reduce noise by using the amplitude of 
this same quantity, namely |ai|, to obtain a less noisy estimate of the power. 

The autocorrelation function at a lag of one is obtained by averaging N-\ complex products, 



The mean Doppler frequency is then commonly obtained from the angle of the resulting phase as 



P 0 = (\/N)Z\z i \ 2 . 



(1) 



tf! = [1/(N-1)]2 z/z l+1 . 



(2) 



/ rf =[PRF/(2a)] arg(a,). 



(3) 
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This is usually displayed at each pixel in blue or red, depending on whether the Doppler shift is 
positive or negative, indicating the direction of flow. 

Defining 

Pi = N (4) 

and 

<h = arg (ai), 

we have 

ai=Pie j ^. 

The invention claimed here displays Pi instead of P a as Power Doppler. The advantage of this 
procedure is that it reduces noise as shown below. 

Consider the case of no signal, i.e., noise only. From equation (1) we see that P 0 results in an 
estimate of the noise variance, = E (|w,| 2 ) while equation (2) reveals that P\ is close to zero if 
the noise is white, i.e., (E (n*n i+ i) = 0). Assuming, at the other extreme, a constant amplitude 
signal with constant Doppler, we have 

zrre j2# + w,-. (5) 

As a first approximation and as the mean or expected value of the actual result, we have 

P 0 * r 2 ifw, = 0 

« o 2 ifr = 0 
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and 

P x * r 2 ifn, = 0 
~ Oifr = 0. 

Since the signal estimate is approximately the same in the absence of noise, Pi is preferred over 
P 0 because it is more immune to noise. 

An alternative to using Pi is to use 

P=[(1/N)SN] 2 . (6) 
Using equation (5), 

P * r 2 ifn, = 0 
« [E|w,|] 2 ifr = 0. 
Let z t = n i = x i + jj^, so that 

\n$\ = <Jxt+yf < 

If x, and y,- are independent normal random variables with zero mean and variance <j\ = /If 
is Rayleigh with mean 

E \n t \ = o x 4itl1 

Hence 

5 



P ~ r 2 if = 0 
~ o 2 7i/4ifr = 0. 

Note that this alternative requires computing N square roots, and yet is almost as noisy as P 0 . 

The advantage of Pi of equations (2) and (4) over P 0 of equation (1) and P of equation (6) can be 
explained in terms of the advantage of coherent signal integration over non-coherent power 
estimation. For a constant signal in white noise, z t = s + n u we estimate the signal as 

1 N 

* = —!>,. (7) 

The signal vectors add coherently while the noise vectors add rms-wise so that E(s) = s 
andE \s\ 2 = s 2 +c 2 / N, where o 2 =E|w/| 2 is the noise power. The above two expected values are 
the means of our estimate of the signal and of our estimate of the signal power, respectively. If 
we were to estimate the signal power by equation (1) instead of by squaring the result of 
equation (7), the mean of our estimate would be s 2 + a 2 ; that is, noise would add up as rapidly 
as signal. 

In Doppler ultrasound, we cannot use equation (7) for coherent signal addition because the signal 
is time varying as, for example, in equation (3), where Sj = r exp(]2nf d iS) where fd is the Doppler 
frequency, / is the discrete time index, and S is the pulse repetition interval (PRI). The signal 
phasor rotates from pulse to pulse so that E(f)= 0. In this instance, the signal adds rms-wise, 
just like the noise. However, if we use equations (2) and (4) we obtain (ignoring cross products 
between signal and noise) 
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E(P,) = E(\a x \) * £(| ^§ r 2 expt/2*/^) + 1) = r 2 + ^ 
where p =0 for white noise. 

This last result comes about because the signal phasors being summed are identical, while the 
random noise phasors add rms-wise. Note that the phase of the constant signal component of a\ 
is 2^PRF, leading to equation (3). 

The disclosure that follows illustrates the exemplary use of Doppler shifted acoustic 
imaging to map and track blood flow. It illustrates the techniques and systems that profit from 
application of the hereindisclosed invention. Alternative probe geometries to which the 
hereindisclosed invention could be profitably applied are illustrated in US Pat. #6,524,253 Bl, 
which is incorporated by reference in its entirety herein. 
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SUMMARY OF THE INVENTION 
There is provided, in accordance with the present invention, a new, useful, and 

10 unobvious method of determining parameters of blood flow, such as vector velocity, 
blood flow volume, and Doppler spectral distribution, using sonic energy (ultrasound) 
and a novel thinned array. Also provided is a novel method of tracking blood flow 
and generating a three dimensional image of a blood vessel of interest that has much 
greater resolution than images produced using heretofore known ultrasound devices 

15 and methods. 

Broadly, the present invention extends to a method for determining a parameter of 
blood flow in a blood vessel of interest, comprising the steps of: 

20 a) providing an array of sonic transducer elements, wherein the element spacing 
in the array is greater than, equal or less than a half wavelength of the sonic energy 
produced by the elements, wherein at least one element transmits sonic energy, and 
a portion of the elements receive sonic energy; 

25 b) directing sonic energy produced by the at least one element of the array into 
a volume of the subject's body having the blood vessel of interest, 

c) receiving echoes of the sonic energy from the volume of the subject's body 
having the blood vessel of interest; 

30 

d) reporting the echoes to a processor programmed to 

i) Doppler process the echoes to determine radial velocity of the blood 
flowing in the blood vessel of interest; 
35 ii) calculate a three dimensional position of blood flow in the vessel of 

interest; and 
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calculate the parameter of blood flow in the blood vessel at the three 
dimensional position calculated in step (ii); and 



(e) displaying the parameter on a display monitor that is electrically connected to 
5 the processor. 

Moreover, a method of the present invention permits an operator examining a subject 
to obtain information on blood flow in a particular region of the blood vessel of 
interest. 

10 

As used herein, the phrases "element spacing" and "distance between the elements" 
can be used interchangeably and refer to the distance between the center of 
elements of an array. 



15 Various methods can be used to determine the three dimensional position of blood 
flow. In a particular embodiment, the method comprises the steps of having the 
processor programmed to: 

i) determine a sum beam, an azimuth difference beam and an elevation 
difference beam from the echoes received from the blood vessel of 

20 interest; 

ii) modulate the directions of the transmitted and received sonic energy 
based upon the sum, azimuth difference and elevation difference 
beams in order to lock on to the highest Doppler energy calculated 
from echoes from the flow of blood in the blood vessel of interest, and 



25 



iii) calculate the three dimensional position of the highest Doppler energy 
from the blood flow in the vessel of interest. 



Optionally, the processor can also be programmed to determine at least one 
30 additional beam having an angle between the azimuth difference beam and the 
elevation difference beam prior to modulating the directions of the transmitted and 
received sonic energy, wherein the at least one additional beam is used to modulate 
the directions of the transmitted and received sonic energy. Naturally, the angle of 
the at least one additional beam can vary. In a particular embodiment, the at least 
35 one additional beam is at an angle that is orthogonal to the blood vessel of interest. 
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Moreover, the present invention extends to a method as described above, wherein 
steps (b) through (e) are periodically repeated so that the three dimensional position 
of blood flow in the vessel of interest is tracked, and the parameter of blood flow is 
periodically calculated and displayed on the display monitor. In a particular 
5 embodiment, the period of time between repeating steps (b) through (e) is sufficiently 
short so that the parameter being measured remains constant, e.g., 20 milliseconds. 

The present invention further extends to a method for determining a parameter of 
blood flow in a particular region of a blood vessel of interest, comprising the steps of: 

10 

a) providing an array of sonic transducer elements, wherein the element spacing 
in the array is greater than, equal or less than a half wavelength of the sonic energy 
produced by the elements, wherein at least one element transmits sonic energy, and 
a portion of the elements receive sonic energy; 

15 

b) directing sonic energy produced by the at least one element of the array into 
a volume of the subjects body having the particular region of the blood vessel of 
interest, 

20 c) receiving echoes of the sonic energy from the volume of the subject's body 
having the particular region of the blood vessel of interest; 

d) reporting the echoes to a processor programmed to 

25 ■) Doppler process the echoes to determine radial velocity of the blood 

flowing in the particular region of the blood vessel of interest; 

ii) calculate a three dimensional position of blood flow in the particular 
region of the blood vessel of interest; and 

iii) calculate the parameter of blood flow in the particular region of the 

30 blood vessel of interest at the three dimensional position calculated in 

step (ii); and 

(e) displaying the parameter on a display monitor that is electrically connected to 
the processor. 



35 
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A particular method of calculating the three dimensional position of blow flow in such 
a method of the present invention comprises having the processor programmed to: 
i) determine a sum beam, an azimuth difference beam and an 
elevation difference beam from the echoes received from the 
particular region of the blood vessel of interest; 

ii) modulate the directions of the transmitted and received sonic energy 
based upon the sum, azimuth difference and elevation difference 
beams in order to lock on to the highest Doppler energy calculated 
from echoes received from the flow of blood in the particular region of 
the blood vessel of interest, and 

iii) calculate the three dimensional position of the highest Doppler energy 
from the blood flow in the particular region of the blood vessel of 
interest. 

As explained above, at least one additional beam can also be determined and used 
to calculate the three dimensional position. 

Furthermore, the present invention extends to a method for determining a parameter 
of blood flow in a blood vessel of interest, comprising the steps of: 

a) providing an array of sonic transducer elements, wherein the element spacing 
in the array is greater than, equal or less than a half wavelength of the sonic energy 
produced by the elements, wherein at least one element transmits sonic energy, and 
a portion of elements receive sonic energy; 

b) directing sonic energy produced by the at least one element of the array into 
a volume of the subject's body having the blood vessel of interest, 

c) receiving echoes of the sonic energy from the volume of the subject's body 
having the blood vessel of interest; 

d) reporting the echoes to a processor electrically connected to the elements of 
the array, wherein the processor is programmed to 
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i) Doppler process the echoes to determine radial velocity of the blood 

flowing in the blood vessel of interest; 

ii) determine a sum beam, an azimuth difference beam and an elevation 
5 difference beam from the echoes received from the blood vessel of 

interest; 

iii) modulate the directions of the transmitted and received sonic energy 
based upon the sum, azimuth difference and elevation difference 

10 beams in order to lock on to the highest Doppler energy calculated 

from echoes from the flow of blood in the blood vessel of interest, 



iv) calculate the three dimensional position of the highest Doppler energy 
from the blood flow in the vessel of interest; and 

15 

v) calculate the parameter of blood flow in the blood vessel at the three 
dimensional position calculated in step (iv); and 



(e) displaying the parameter on a display monitor that is electrically connected to 
20 the processor. 



As explained above, an operator performing a method of the present invention can 
obtain blood flow parameters from a blood vessel of interest, and even from a 
particular region of a blood vessel of interest. 

25 

Moreover, the present invention extends to a method for determining a parameter of 
blood flow in a particular region of a blood vessel of interest, comprising the steps of: 

a) providing an array of sonic transducer elements, wherein the element spacing 
30 in the array is greater than, equal or less than a half wavelength of the sonic energy 

produced by the elements, wherein at least one element transmits sonic energy, and 
a portion of the elements receive sonic energy; 

b) directing sonic energy produced by the at least one element of the array into 
35 a volume of the subject's body having the particular region of the blood vessel of 

interest, 
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c) receiving echoes of the sonic energy from the volume of the subject's body 
having the particular region of blood vessel of interest; 

5 d) reporting the echoes to a processor electrically connected to the elements of 
the array, wherein the processor is programmed to 

i) Doppler process the echoes to determine radial velocity of the blood 
flowing in the particular region of the blood vessel of interest; 

10 

ii) determine a sum beam, an azimuth difference beam and an elevation 
difference beam from the echoes received from the particular region of 
the blood vessel of interest; 

15 iii) modulate the directions of the transmitted and received sonic energy 

based upon the sum, azimuth difference and elevation difference 
beams in order to lock on to the highest Doppler energy calculated 
from echoes from the flow of blood in the particular region of the blood 
vessel of interest, 

20 

iv) calculate the three dimensional position of the highest Doppler energy 
from the blood flow in the particular region of the blood vessel of 
interest; and 

25 v) calculate the parameter of blood flow in the particular region of the 

blood vessel at the three dimensional position calculated in step (iv); 
and 



(e) displaying the parameter on a display monitor that is electrically connected to 
30 the processor. 

In another embodiment, the present invention extends to a device for determining a 
parameter of blood flow in a blood vessel of interest, comprising: 



35 



a) an array of sonic transducer elements, wherein the element spacing in the 
array is greater than, equal or less than a half wavelength of the sonic energy 
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produced by the elements, and at least one element transmits sonic energy, and a 
portion of the elements receive sonic energy; 

b) a processor electrically connected to the array so that echoes received from a 
5 volume of the subject's body having the blood vessel of interest due to directing sonic 
energy produced by the at least one element of the array into the subject's body is 
reported to the processor, wherein the processor is programmed to: 

i) Doppler process the echoes to determine radial velocity of the blood 
flowing in the blood vessel of interest; 
10 ii) calculate a three dimensional position of blood flow in the blood vessel 

of interest; and 

iii) calculate the parameter of blood flow in the blood vessel of interest at 
the three dimensional position calculated in step (ii); and 

15 (c) a display monitor that is electrically connected to the processor which 
displays the parameter of blood flow calculated by the processor. 

A parameter of blood that can be determined with a device of the present invention 
includes blood flow volume, vector velocity, Doppler spectral distribution, etc. The 
20 parameter being measured can be an instantaneous value, or an average value 
determined over a heart cycle. 

Moreover, the present invention extends to a device as described above, wherein 
the processor is programmed to: 

25 

i) 



30 ii) 



determine a sum beam, an azimuth difference beam and an elevation 
difference beam from the echoes received from the blood vessel of 
interest after Doppler processing the echoes; 

modulate the directions of the transmitted and received sonic energy 
based upon the sum, azimuth difference and elevation difference 
beams in order to lock on to the highest Doppler energy calculated 
from echoes from the flow of blood in the blood vessel of interest, 

iii) calculate the three dimensional position of the highest Doppler 
energy from the blood flow in the vessel of interest; and 
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iv) calculate the parameter of blood flow in the blood vessel of interest at 

the three dimensional position calculated in (iii). 

Optionally, a processor of a device of the present invention can be further 
5 programmed to determine at least one additional beam having an angle between the 
azimuth difference beam and the elevation difference beam prior to modulating the 
directions of the transmitted and received sonic energy, wherein the at least one 
additional beam is used to modulate the directions of the transmitted and received 
sonic energy. In a particular embodiment, the at least one additional beam is at an 
10 angle that is orthogonal to the blood vessel of interest. 

Moreover, in a another embodiment of a device of the present invention, the distance 
between the elements of the array is greater than 14 the wavelength of the sonic 
energy generated by the at least one element. 

15 

Furthermore, the present invention extends to a device for determining a parameter 
of blood flow in a blood vessel of interest, comprising: 

a) an array of sonic transducer elements, wherein the element spacing in the 
20 array is greater than, equal or less than a half wavelength of the sonic energy 

produced by the elements, and at least one element transmits sonic energy, and 
portion of the elements receive sonic energy; 

b) a processor electrically connected to the array so that echoes received from a 
25 volume of the subject's body having the blood vessel of interest due to directing sonic 

energy produced by the at least one element of the array into the subject's body is 
reported to the processor, wherein the processor is programmed to: 

i) Doppler process the echoes to determine radial velocity of the blood 
flowing in the blood vessel of interest; 
30 ii) calculate a three dimensional position of blood flow in the blood vessel 

of interest; and 

iii) calculate the parameter of blood flow in the blood vessel of interest at 
the three dimensional position calculated in step (ii) 

35 (c) a display monitor that is electrically connected to the processor which 
displays the parameter of blood flow calculated by the processor. 
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Particular parameters of blood flow that can be determined with a device of the 
present invention include, but certainly are not limited to blood flow volume, vector 
velocity, and Doppler spectral distribution. The parameter being measured can be an 
5 instantaneous value, or an average value determined over a heart cycle. 

In addition, a processor of a device of the present invention can be further 
programmed to determine at least one additional beam having an angle between the 
azimuth difference beam and the elevation difference beam prior to modulating the 
10 directions of the transmitted and received sonic energy, wherein the at least one 
additional beam is used to modulate the directions of the transmitted and received 
sonic energy. In a particular embodiment, the at least one additional beam is at an 
angle that is orthogonal to the blood vessel of interest. 

15 Moreover, the present invention extends to. a method for generating a three 

dimensional image using sonic energy of a blood vessel of interest in a subject, the 
method comprising the steps of: 

a) providing an array of sonic transducer elements, wherein the element spacing 
20 in the array is greater than, equal or less than a half wavelength of the sonic energy 

produced by the elements, wherein at least one element transmits sonic energy, and 
a portion of the elements receive sonic energy; 

b) directing sonic energy produced by the at least one element of the array into 
25 a volume of the subject's body having the blood vessel of interest, 

c) receiving echoes of the sonic energy from the volume of the subject's body 
having the blood vessel of interest; 

30 d) reporting the echoes to a processor programmed to 

i) Doppler process the echoes to determine radial velocity of the blood 
flowing in the blood vessel of interest; 

ii) calculate a three dimensional position of blood flow in the blood vessel 
35 of interest; 
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iii) repeat steps (i) through (ii) to generate a plurality of calculated three 

dimensional positions; and 
vi) generate a three dimensional image of the blood vessel of interest 

from the plurality of calculated three dimensional positions; and 

5 

(e) displaying the three dimensional image on a display monitor that is electrically 
connected to the processor. 

Furthermore, the present invention permits an operator utilizing a method of the 
10 present invention to generate a three dimensional image of not only a blood vessel in 
the body, but even a particular region of a blood vessel in the body. 

Numerous means available for calculating the three dimensional position of a blood 
vessel and even a particular portion of a blood vessel are encompassed by the 
15 present invention. A particular means comprises having the programmed processor: 

i) determine a sum beam, an azimuth difference beam and an elevation 
difference beam from the echoes received from the blood vessel of 
interest after Doppler processing the echoes; 



20 



25 



ii) modulate the directions of the transmitted and received sonic energy 
based upon the sum, azimuth difference and elevation difference 
beams in order to lock on to the highest Doppler energy calculated 
from echoes from the flow of blood in the blood vessel of interest, and 

iii) calculate the three dimensional position of the highest Doppler energy 
from the blood flow in the vessel of interest, and 



iv) repeat steps (i) through (iii) to generate a plurality of calculated three 
30 dimensional positions. 

Optionally, a processor of a method of the present invention can also be 
programmed to determine at least one additional beam having an angle between the 
azimuth difference beam and the elevation difference beam prior to modulating the 
35 directions of the transmitted and received sonic energy, and the at least one 
additional beam is also used to modulate the directions of the transmitted and 
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received sonic energy, and calculate the three dimensional position of the highest 
Doppler energy. In a particular embodiment, the at least one additional beam is at 
an angle that is orthogonal to the blood vessel of interest. 



5 The present invention also extends to a method for generating a three dimensional 
image of a blood vessel of interest in a subject using sonic energy, the method 
comprising the steps of: 



a) providing an array of sonic transducer elements, wherein the element spacing 
10 in the array is greater than, equal or less than a half wavelength of the sonic energy 
produced by the elements, wherein at least one element transmits sonic energy, and 
a portion of the elements receive sonic energy; 



b) directing sonic energy produced by the at least one element of the array into 
15 a volume of the subject's body having the blood vessel of interest, 

c) receiving echoes of the sonic energy from the volume of the subject's body 
having the blood vessel of interest; 



20 d) reporting the echoes to a processor programmed to 



i) Doppler process the echoes to determine radial velocity of the blood 
flowing in the blood vessel of interest; 



25 ii) determine a sum beam, an azimuth difference beam and an elevation 

difference beam from the echoes received from a portion of the blood 
vessel of interest; 

iii) modulate the directions of the transmitted and received sonic energy 
30 based upon the sum, azimuth difference and elevation difference 

beams in order to lock on to the highest Doppler energy calculated 
from echoes from the flow of blood in the blood vessel of interest, 



calculate the three dimensional position of the highest Doppler energy 
from the blood flow in the vessel of interest; and 
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v) repeat steps (i) through (iv) to generate a plurality of calculated three 

dimensional positions; 

vi) generate a three dimensional image of the blood vessel of interest 
5 from the plurality of calculated three dimensional positions; and 

(e) displaying the three dimensional image on a display monitor that is electrically 
connected to the processor. 

10 Optionally, the three dimensional image can be of a particular region of a blood 
vessel of interest. Moreover, a processor of a method described herein can also 
determine at least one additional beam having an angle between the azimuth 
difference beam and the elevation difference beam prior to modulating the directions 
of the transmitted and received sonic energy, and the at least one additional beam is 

15 also used to modulate the directions of the transmitted and received sonic energy, 
and calculate the three dimensional position of the highest Doppler energy. Angles 
for use with the at least one additional beam are described above. 

Moreover, in another embodiment of the present invention, the distance between the 
20 elements of the array is greater than Vi the wavelength of the sonic energy generated 
by the at least one element. 

Furthermore, the present invention extends to a device generating a three 
dimensional image of a blood vessel of interest in a subject using sonic energy, 
25 comprising: 

a) an array of sonic transducer elements, wherein the element spacing in the 
array is greater than, equal or less than a half wavelength of the sonic energy 
produced by the elements, and at least one element transmits sonic energy, and a 

30 portion of the elements receive sonic energy; 

b) a processor electrically connected to the array so that echoes received from a 
volume of the subject's body having the blood vessel of interest due to directing sonic 
energy produced by the at least one element of the array into the subject's body is 

35 reported to the processor, wherein the processor is programmed to: 
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i) Doppler process the echoes to determine radial velocity of the blood 

flowing in the blood vessel of interest; 
ii) calculate a three dimensional position of blood flow in the blood vessel 
of interest; 

5 iii) repeat steps (i) through (ii) to generate a plurality of calculated three 

dimensional positions; 
v) generate a three dimensional image from the plurality of calculated 
three dimensional positions, and 

10 (c) a display monitor that is electrically connected to the processor which 
displays the three dimensional image. 

As explained above, a device of the present invention permits an operator to 
generate and display three dimensional images of a blood vessel of interest, and 
15 even of a particular region of a blood vessel that the operator wants to investigate 
closely. Moreover, in a particular embodiment, a processor of a device of the present 
invention can be programmed to calculate the three dimensional position of a blood 
vessel by 

20 i) determining a sum beam, an azimuth difference beam and an 

elevation difference beam from the echoes received from the blood 
vessel of interest after Doppler processing the echoes; 

ii) modulating the directions of the transmitted and received sonic energy 
25 based upon the sum, azimuth difference and elevation difference 

beams in order to lock on to the highest Doppler energy calculated 
from echoes from the flow of blood in the blood vessel of interest, 

iii) calculating the three dimensional position of the highest Doppler 
30 energy from the blood flow in the vessel of interest; and 

iv) repeat steps (I) through (iii) in order to generate a plurality of 
calculated three dimensional positions used to generate the three 
dimensional image. 

35 
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Optionally, the processor can be programmed to further determine at least one 
additional beam having an angle between the azimuth difference beam and the 
elevation difference beam prior to modulating the directions of the transmitted and 
received sonic energy, wherein the at least one additional beam is used to modulate 
5 the directions of the transmitted and received sonic energy. The angle between the 
azimuth difference beam and the elevation difference beam of the additional beam 
can vary. In a particular embodiment, the at least one additional beam is at an angle 
that is orthogonal to the blood vessel of interest. 

10 Furthermore, the present invention extends to a thinned array for use in an 

ultrasound device, comprising a plurality of sonic transducer elements, wherein the 
element spacing in the array is greater than a half wavelength of the sonic energy 
produced by the elements, and the elements are positioned and sized within the 
array, and sonic energy is electronically steered by the elements so that any grating 

15 lobes produced by the sonic energy are suppressed. In a particular embodiment, 
the elements positioned and sized so that they are flush against each other. 

Hence, the current invention performs blood velocity monitoring by collecting Doppler 
data in three dimensions; azimuth, elevation, and range (depth); so that the point (in 
20 three dimensional space) at which the velocity is to be monitored can be acquired 
and tracked when the patient or the sensor moves. The invention also produces a 
three dimensional map of the blood flow and converts measured radial velocity to 
true vector velocity. 

25 Moreover, in this invention, once the desired signal is found, it will be precisely 
located and continually tracked with accuracy far better than the resolution. A 
heretofore unknown method to achieve sub-resolution tracking and mapping involves 
a novel and unobvious extension of a procedure called "monopulse". Monopulse 
tracking has been used in military applications for precisely locating and tracking a 

30 point target with electromagnetic radiation. However, it has never been utilized in 
connection with sonic waves to determine the velocity of moving fluids in vivo. 



35 



This invention provides: (1) affordable three-dimensional imaging of blood flow using 
a low-profile easily-attached transducer pad, (2) real-time vector velocity, and (3) 
long-term unattended Doppler-ultrasound monitoring in spite of motion of the patient 
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or pad. None of these three features are possible with current ultrasound 
equipment or technology. 

The pad and associated processor collects and Doppler processes ultrasound blood 
velocity data in a three-dimensional region through the use of a two-dimensional 

5 phased array of piezoelectric elements on a planar, cylindrical, or spherical surface. 
Through use of unique beamforming and tracking techniques, the invention locks 
onto and tracks the points in three-dimensional space that produce the locally 
maximum blood velocity signals. The integrated coordinates of points acquired by 
the accurate tracking process is used to form a three-dimensional map of blood 

10 vessels and provide a display that can be used to select multiple points of interest for 
expanded data collection and for long term continuous and unattended blood flow 
monitoring. The three dimensional map allows for the calculation of vector velocity 
from measured radial Doppler. 

15 In a particular embodiment, a thinned array (greater than half-wavelength element 
spacing of the transducer array) is used to make a device of the present invention 
inexpensive and allow the pad to have a low profile (fewer connecting cables for a 
given spatial resolution). The array is thinned without reducing the receiver area by 
limiting the angular field of view. The special 2-D phased array used in this invention 

20 makes blood velocity monitoring inexpensive and practical by (1) forming the beams 
needed for tracking and for re-acquiring the blood velocity signal and by (2) allowing 
for an element placement that is significantly coarser than normal half-wavelength 
element spacing. The limited range of angles that the array must search allows for 
much less than the normal half wavelength spacing without reducing the total 

25 receiver area. 

Grating lobes due to array thinning can be reduced by using wide bandwidth and 
time delay steering. The array, or at least one element of the array, is used to 
sequentially insonate the beam positions. Once the region of interest has been 

30 imaged and coarsely mapped, the array is focused at a particular location on a 
particular blood vessel for measurement and tracking. Selection of the point or 
points to be measured and tracked can be based on information obtained via 
mapping and may be user guided or fully automatic. Selection can be based, for 
example, on peak response within a range of Doppler frequencies at or near an 

35 approximate location. 
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In the tracking mode a few receiver beams are formed at a time: sum, azimuth 
difference, elevation difference, and perhaps, additional difference beams, at angles 
other than azimuth ( = 0 degrees) and elevation (= 90 degrees). Monopulse is 
applied at angles other than 0 and 90 degrees (for example 0, 45, 90, and 135 
5 degrees) in order to locate a vessel in a direction perpendicular to the vessel. When 
the desired (i.e. peak) blood velocity signal is not in the output, this is instantly 
recognized (e.g., a monopulse ratio, formed after Doppler filtering, becomes non- 
zero) and the array is used to track (slow movement) or re-acquire (fast movement) 
the desired signal. Re-acquisition is achieved by returning to step one to form and 

10 Doppler-process a plurality of beams in order to select the beam (and the time delay 
or "range gate") with the most high-Doppler (high blood velocity) energy. This is 
followed by post-Doppler monopulse tracking to lock a beam and range gate on to 
the exact location of the peak velocity signal. In applications such as transcranial 
Doppler, where angular resolution based on wavelength and aperture size is 

15 inadequate, fine mapping is achieved, for example, by post-Doppler monopulse 
tracking each range cell of each vessel, and recording the coordinates and 
monopulse-pair angle describing the location and orientation of the monopulse null. 
With a three-dimensional map available, true vector velocity can be computed. For 
accurate vector flow measurement, the monopulse difference is computed in a 

20 direction orthogonal to the vessel by digitally rotating until a line in the azimuth- 
elevation or C-scan display is parallel to the vessel being monitored. The aperture is 
more easily rotated in software (as opposed to physically rotating the transducer 
array) if the aperture is approximately circular (or eliptical) rather than square (or 
rectangular). Also, lower sidelobes result by removing elements from the four 

25 corners of a square or rectangular array in order to make the array an octagon. 

In this invention, as long as (1) a blood vessel or (2) a flow region of a given velocity 
can be resolved by finding a 3-D resolution cell through which only a single vessel 
passes, that vessel or flow component can then be very accurately located within the 
cell. Monopulse is merely an example of one way to attain such sub-resolution 
30 accuracy (SRA). Other methods involve "super-resolution" or "parametric'' 

techniques used in "modern spectral estimation", including the MUSIC algorithm and 
autoregressive modeling, for example. SRA allows an extremely accurate map of 3- 
D flow. 

Furthermore, the present invention utilizes post-Doppler, sub-resolution tracking and 
35 mapping; it does Doppler processing first and uses only high Doppler-frequency data. 
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This results in extended targets since the active vessels approximate "lines" as 
opposed to "points". In three-dimensional space, these vessels are resolved, one 
from another. At a particular range, the monopulse angle axis can be rotated (in the 
azimuth-elevation plane) so that the "line" becomes a "point" in the monopulse angle 
5 direction. That point can then be located by using super-resolution techniques or by 
using a simple technique such as monopulse. By making many such measurements 
an accurate 3-D map of the blood vessels results. 

Methods for extending the angular field of view of the thinned array (that is limited by 
grating lobes) include (1) using multiple panels of transducers with multiplexed 
10 processing channels, (2) convex V-shaped transducer panels, (3) cylindrical shaped 
transducer panel, (4) spherical shaped transducer panel, and (5) negative ultrasound 
lens. If needed, moving the probe and correlating the sub-images can create a map 
of an even larger region. 

Active digital beamforming can also be utilized, but the implementation depends on a 
15 choice to be made between wideband and narrowband implementations. If 

emphasis is on high resolution mapping of the blood vessels, then a wide bandwidth 
(e.g., 50% of the nominal frequency) is used for fine range resolution. If emphasis is 
on Doppler spectral analysis, measurement, and monitoring, the map is only a tool. 
In this case, a narrowband, low cost, low range-resolution, high sensitivity 
20 implementation might be preferred. A wideband implementation would benefit in 
performance (higher resolution, wider field of view, and reduced grating lobes) using 
time-delay steering while a narrowband implementation would benefit in cost using 
phase-shift steering. The invention can thus be described in terms of two preferred 
implementations. 

25 In a wideband implementation, time delay steering can be implemented digitally for 
both transmit and receive by over-sampling and digitally delaying in discrete sample 
intervals. In a narrowband implementation, (1) phase steering can be implemented 
digitally (digital beamforming) for both transmit and receive, and (2) bandpass 
sampling (sampling at a rate lower than the signal frequency) can be employed with 

30 digital down-conversion and filtering. 

Accordingly, it is an object of the present invention to locate the point in three 
dimensional space having the greatest high-Doppler energy, and determining 
coordinates for that point. With that information, and the radial velocity of the blood 
35 flowing through the blood vessel at that point, a variety of blood flow parameters can 
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be calculated at that point, including, but not limited to vector velocity of blood flow, 
volume of blood flow, or Doppler spectral distribution. The parameter being 
measured can be an instantaneous value, or an average value determined over a 
heart cycle. 

5 

It is also an object of the present invention to continuously track and map in vivo the 
point in three dimensional space having the greatest Doppler-energy, and using the 
coordinates to generate a three dimensional image of a blood vessel and blood flow 
therein that possess a much greater resolution than images generated using 
10 heretofore known Doppler ultrasound methods and devices. 

It is yet another object of the present invention to provide a thinned array which does 
not utilize the number of element transducers as are required with heretofore known 
Doppler ultrasound devices. As a result, the decreased number of elements in the 
15 array decreases size of the array utilized and provides a patient being analyzed with 
mobility that would not be available if using conventional ultrasound devices to obtain 
blood flow parameters such as vector velocity, blood flow volume, and Doppler 
spectral distribution. The parameter being measured can be an instantaneous value, 
or an average value determined over a heart cycle. 

20 

These and other aspects of the present invention will be better appreciated by 
reference to the following drawings and Detailed Description. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates the Blood Flow Mapping Monitor in use with a Transcranial Doppler 
25 Probe, as an example. 

FIG.2 shows a 64-element bistatic ultrasound transducer array example, where, with 
D = 2d, the same elements are reconfigured differently for transmit and receive 
during the acquisition phase of operation. Fig. 2 (a) shows the Receive 
Configuration, where all 64 elements receive at once. Fig. 2 (b) shows the Transmit 
30 Configuration, where, during acquisition, the 16 sub-apertures transmit one at a time. 

FIG. 3 is an example overall block diagram of a blood flow mapping monitor 
embodiment. 
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FIG. 4 illustrates ultrasound beam coverage for the TCD array example of Fig. 2. 
The left illustration shows 25 digitally beam-formed beams, as an example. On the 
right, is shown, for that example, the manner in which the transmit beam 
encompasses 21 receive beams in the acquisition mode. 

5 FIG. 5 shows one-dimensional patterns for a bistatic transducer array with D = 2 d as 
in FIG. 2. FIG. 5a (top) shows the transmit element pattern. FIG. 5b shows the 
receive Element Pattern and Array Pattern with the receiver beam steered to 
broadside (x=0). The Array Pattern has Grating Lobes (Receiver Ambiguities). FIG. 
5c shows the resultant two-way beam pattern (product of all three patterns above). 

10 The Grating Lobes are suppressed. 

FIG. 6 is the same as FIG. 5, with the receive array beam steered to x = 0.2. 

FIG. 7 shows the Two-way pattern of a receiver beam steered to the half power point 
(x=0.2). This is Fig. 6c plotted in dB. 

FIG.8 shows a one-dimensional representation of the example of FIG. 4. FIG. 8a 
15 shows the product of transmit and receive Element Patterns. FIG. 8b plots a set of 
five receive beams showing Grating Lobes of the Thinned Array. FIG. 8c plots the 
resultant two-way beams with Grating Lobes suppressed. 

FIG.9 is a block diagram of one possible embodiment of the Transmit-Receive 
Electronics for a Bistatic Ultrasound Imaging Sensor and Blood Monitoring Monitor. 

20 FIG. 10 shows the receiver channel signal spectrum illustrating functions performed 
by the FPGA of Fig. 9 on each of the 64 received signals for a narowband case. 

FIG.1 1 shows the geometry involved in using azimuth monopulse to more accurately 
determine the cross-range location of a vessel. The range resolution is better than 
the cross-range resolution and the measured radial velocity field or color flow map 
25 has been utilized to rotate and orient the azimuth and elevation axes so that the 
center of the vessel is vertical, at approximately zero azimuth. The black circular 
cylinder represents the location of all points within the spatial resolution cell that have 
a particular velocity. 



30 FIG. 12 shows the geometry involved in using Doppler ultrasound to determine the 
diameter of a vessel or the velocity field within the vessel. While the initial 3-D 
orientation of the vessel is general, a measured 3-D radial velocity field or 3-D color 
flow map has been utilized to rotate and orient the azimuth and elevation axes so 
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that the center of the vessel is vertical, at approximately zero azimuth. In other 
words, the coordinate system has been rotated about the depth-axis so that the 
centerline of the vessel is in the depth-elevation plane. This can be accomplished 
either by a change of coordinates in software or by physically rotating the ultrasound 
5 probe. The black circular cylinder represents the location of all points within the 
illustrated box that have a particular velocity. The diameter of the cylinder is then 
measured as the azimuth extent of a high-resolution depth-azimuth or B-scan image 
at the Doppler frequency under examination. 

FIG. 13 illustrates the Blood Flow Mapping Monitor in use with a Transcranial 
10 Doppler Probe, as an example. 

FIG. 14 shows a 52-element ultrasound transducer array example, based on an 8 by 
8 rectangular array of elements with 3 elements removed from each corner to make 
the array octagonal instead of rectangular or square. For this example, the elements 
are square (di = d 2 = d) and Ud = 8. 

15 FIG. 15 shows a typical pattern of electronically scanned beams produced by the 
array in FIG. 14. The beam width is nominally, given by the signal wavelength 
divided by the size, L, of the array. The angular field of view (F.O.V.) is limited by the 
maximum angle to which the array can be steered without producing grating lobes 
that are not sufficiently attenuated by the pattern of the individual dxd element. 

20 FIG. 16 shows one-dimensional patterns for an eight-element monostatic linear 
transducer array corresponding to a column or a row in FIG. 16. FIG. 16a (top) 
shows the Element Pattern and Array Pattern with the beam steered to broadside 
(x=0). The Array Pattern has Grating Lobes (Receiver Ambiguities). FIG. 16b shows 
the resultant beam pattern. The Grating Lobes are suppressed. 

25 FIG. 17 is the same as FIG. 16, with the array beam steered to an angle at which a 
grating lobe exceeds the highest sidelobe. The thinned array of Figure 16 should not 
be steered beyond ± arcsin (X/5d) (± 4.7° for the example used) if grating lobes are to 
be suppressed. 

FIG. 18 shows the pattern of a beam steered to the point where the grating lobe 
30 problem appears. This is Fig. 17b plotted in dB. 

FIG. 19 shows a dual 52-active-element ultrasound transducer array example (similar 
to that in FIG. 14) with a total of 1 16 elements, 52 of which are used at a time. 
Figure 19 B shows that the two sub-arrays are in two different planes, tilted to reduce 
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the overlap between beams from the two sub-arrays and maximize the azimuth 
angular field of view. 

FIG. 20 shows a 52-active-element ultrasound transducer array example (similar to 
that in FIG. 14) with a total of 84 elements (52 of which are used at a time) and with a 
5 slightly convex cylindrical shape. The indicated L, x L 2 ' sub-aperture would be 
activated for the formation of beams pointed to one side. 

FIG. 21 is an example overall block diagram of a blood flow mapping monitor 
embodiment. 

FIG. 22 is a block diagram of one possible embodiment of the analog Transmit- 
10 Receive Electronics for an Ultrasound Imaging Sensor and Blood Monitor. 

FIG. 23 shows the geometry involved in using azimuth monopulse to more 
accurately determine the cross-range location of a vessel. The measured radial 
velocity field or color flow map has been utilized to rotate and orient the azimuth and 
elevation axes so that the center of the vessel is vertical, at approximately zero 
15 azimuth. The black circular cylinder represents the location of all points within the 
spatial resolution cell that have a particular velocity. 

DETAILED DESCRIPTION OF THE INVENTION 

The invention involves (1 ) a family of ultrasound sensors, (2) the interplay of a set of 
core technologies that are unique by themselves, and (3) a number of design options 
20 which represent different ways to implement the invention. To facilitate an 

organizational understanding of this many-faceted invention, a discussion of each of 
the three topics above follows. 

The sensors addressed are all two-dimensional (i.e., planar or on the surface of a 
convex shape such as a section of a cylinder) arrays of piezoelectric crystals for use 

25 in active, non-invasive, instantaneous (or real-time), three-dimensional imaging and 
monitoring of blood flow. The sensors use a unique approach to 3-D imaging of 
blood velocity and blood flow that (1) allows for finer image resolution than would 
otherwise be possible with the same hardware complexity (number of input cables 
and associated electronics) and (2) allows for finer accuracy than would ordinarily be 

30 possible based on the resolution. The invention measures and monitors 3-D vector 
velocity rather than merely the radial component of velocity. 

Moreover, the present invention also utilizes (1) array thinning with large elements 
and limited scanning, (2) array shapes to reduce peak sidelobes and extend the field 
of coverage, (3) post-Doppler sub-resolution tracking, (4) post-Doppler sub-resolution 
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mapping, (5) additional methods for maximizing the angular field of view, and (6) 
various digital beamforming procedures for implementing the mapping, tracking, and 
measurement processes. The present invention also extends to array thinning, 
where the separation between array elements is significantly larger than half the 

5 wavelength. This reduces the number of input cables and input signals to be 

processed while maintaining high resolution and sensitivity and avoiding ambiguities. 
In a transcranial Doppler application, for example, where signal to noise and hence 
receiver array area is of paramount importance, array thinning is possible without 
reducing the receiver array area because a relatively small (compared to other 

10 applications) angular field of view is needed. 

Thinning with full aperture area imposes limitations on the angular field of view. 
Methods for expanding the field of view include using more elements than are active 
at any one time. For example, if the electronics are switched between two identical 
panels, the cross-range field of view at any depth is increased by the size of the 
15 panel. If the panels are pointed in slightly different directions so that overlapping or 
redundant beams are avoided, the field of view is doubled. A generalization of this 
approach involves the use of an array on a cylindrical or spherical surface. 

Once a section of a blood vessel is resolved from other vessels in Doppler, depth, 
and two angles (az and el), Post-Doppler sub-resolution processing locates that 
20 section to an accuracy that is one-tenth to one-twentieth of the resolution. This 
allows for precise tracking and accurate mapping. Tracking provides for the 
possibility of unattended long term monitoring and mapping aids the operator in 
selecting the point or points to be monitored. 

Furthermore, methods of the present invention permit non-invasive, continuous, 
25 unattended, volumetric, blood vessel tracking, ultrasound monitoring and diagnostic 
device for blood flow. It will enable unattended and continuous blood velocity 
measurement and monitoring as well as 3-dimensional vascular tracking and 
mapping using an easily attached, electronically steered, transducer probe that can 
be in the form of a small pad for monitoring application, when desired. Moreover, a 
30 device and method of the present invention have applications in measuring the 
parameters described above in any part of the body. A nonlimiting example 
described below involves a cranial application. However as set forth, a device and 
method of the present have applications in any part of the body, and can be used to 
track and map any blood vessel in the body. A device of the present invention can, 
35 for example: 
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1. Measure and continuously monitor blood velocity with a small low-profile probe 
that can be adhered, lightly taped, strapped, banded, or otherwise easily 
attached to the portion of the body where the vascular diagnosis or monitoring 
is required. 

2. Track and maintain focus on multiple desired blood vessels in spite of 
movement. 

3. Map 3-D blood flow; e.g., in the Circle of Willis (the central network of arteries 
that feeds the brain) or other critical vessels in the cranial volume. 

4. Perform color velocity imaging and display a 3-D image of blood flow that is 
rotated via track ball or joystick until a desired view is selected. 

5. Form and display a choice of projection, slice, or perspective views, including 
(1) a projection on a depth-azimuth plane, a B-scan, or a downward -looking 
perspective, (2) a projection on an azimuth-elevation plane, a C-scan, or a 
forward-looking perspective, or (3) a projection on an arbitrary plane, an 
arbitrary slice, or an arbitrary perspective. 

6. Use a track ball and buttons to position circle markers on the points were 
measurement or monitoring of vector velocity is desired. 

7. Move the track location along the blood vessel by using the track ball to slide 
the circle marker along the image of the vessel. 

8. Display actual instantaneous and/or average vector velocity and/or estimated 
average volume flow. 

9. Maintain a multi-day history and display average blood velocity versus time 
for each monitored vessel over many hours. 

10. Sound an alarm when maximum or minimum velocity is exceeded or when 
emboli count is high; and maintain a log of emboli detected. 

1 1 . Track, map, and monitor small vessels (e.g., 1mm in diameter), resolve 
vessels as close as 4 mm apart (for example), and locate them with an 
accuracy of ±0.1 mm, for example. 

Moreover, as explained herein, numerous methods have applications in obtaining the 
three dimensional coordinates of points along a blood vessel from echoes returned 
from the body, and are encompassed by the present invention. A particular 
nonlimiting example of such a method having applications herein is a novel and 
unobvious variation of monopulse tracking. For tracking purposes utilizing 
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monopulse, up to nine beams are simultaneously formed for each transmit beam 
position. In addition to the "sum" beam that corresponds to the transmitted beam, 
there will either be 4 monopulse difference beams or there will be 8 overlapping 
focused beams. If a cluster of eight focused beams is used, these will be highly 

5 overlapped with the sum beam, and displaced very slightly from the sum beam, with 
their centers equally spaced on a small circle around the center of the sum beam. 
These satellite beams would then operate in pairs to form four difference beams. For 
example, the azimuth Monopulse ratio can be produced in two different ways, which 
will call "liner" and "non-linear". The non-linear method will determine the 

10 magnitudes or the powers of three received signals, left, right, and sum (L, R, and S), 
and compute M a = (|L| - |R|) / |S|. The linear method uses complex signals and 
computes the azimuth monopulse ratio as the real part of the ratio D a / S, where D a = 
L - R. D a is the azimuth difference. 

15 For an ideal point target, the linear method for computing M a results in an excellent 
estimate of the azimuth angle error. It also has the advantage of only requiring 4, 
instead of 8 auxiliary beams. These 4 beams would be an azimuth difference beam, 
D a , an elevation difference beam, and two diagonal difference beams. The individual 
beams, such as L and R, are not needed. However, beam shapes will be highly 

20 distorted by refraction through bone and tissue, and a "sub-optimum" non-linear 
approach might be more robust. 

Regardless of which monopulse method is used, the conventional two difference 
beams used in radar (azimuth difference and elevation difference) may not be 

25 enough. The projection of the high-velocity data on a plane perpendicular to the 
transducer line of sight (the C-scan) will usually be a line, not a point. With multiple 
difference beams, equally spaced in angle, one will be approximately perpendicular 
to the C-scan projection of the vessel. The system will select the monopulse 
difference output with the largest magnitude. This provides an approximate 

30 orientation of the C-scan projection of the vessel. The corresponding monopulse 
ratio (provided the sum beam power exceeds a threshold) is used to correctly re- 
steer and maintain a beam precisely centered on that vessel. 

If the power map output of a Wall filter is used for the monopulse beams, the beam 
35 outputs are power and hence a complex ratio is not available. In that case the 
nonlinear method would be used. An alternative is to use the complex wall filter 
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10 



15 



20 



output, before computing the power, with the linear method. During measurement, 
however, the output of a particular (high velocity) FFT Doppler bin may be used for 
monopulse (provided that the magnitude or power of he sum beam at that Doppler 
exceeds a threshold). In that case either the linear or the nonlinear monopulse ratio 
may be used. 

Another alternative is to use FFT processing and form the monopulse ratio (linearly 
or non-linearly) at the output of a high-velocity Doppler-frequency cell with high sum- 
beam power. For example, set a power threshold and select the highest (positive or 
negative) velocity cell with power that exceeds the threshold. Since the data in a 
single FFT cell is expected to be noisy, this procedure is recommended for a 
measurement dwell, where enough time is spent in a single beam position to have 
both useable velocity resolution and the ability to make several measurements 
(multiple FFTs per frame). 

FFT-Based Monopulse and Monopulse Averaging During Measurement 
In a K pulse dwell, let K=K 1 xK 2l where K y is the number of input pulses used in the 
FFT and K 2 is the number of FFTs. Instead of performing monopulse to re-steer the 
beam every K t pulses, we compute the monopulse ratio at the output of a desired 
high velocity Doppler bin, and average its value over K 2 FFTs. This reduces the 
steering noise while assuring that we are locating the center of the vessel (the 
highest Doppler Energy). We chose the highest Doppler frequency for which the 
minimum sum beam power exceeds a threshold, and utilize only that Doppler cell for 
monopulse. The average is best performed as a weighted average. For example, if 
D n and S n are (say, elevation) difference-beam and sum beam outputs in the nth FFT 
for the selected Doppler bin, we chose: 



M = 




where M n 

n 



K*{D H /S H } 



or M n 



\s} 
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depending on whether linear or non-linear monopulse is used. For linear 
monopulse it might be best to use only one large FFT (K 2 =1 ). For non-linear 
monopulse, the expression simplifies to: 

K 2 

2>J 2 

M = *i? 

[Note that because a ratio is involved (so that beam pointing error is not confused 
with signal strength) even the "linear" method is non-linear.] 

A device of the present invention will allow a person with little training to apply the 
10 sensor and position it based on an easily understood ultrasound image display. The 
unique sensor can continuously monitor artery blood velocity and volume flow for 
early detection of critical events. It will have an extremely low profile for easy 
attachment, and can track selected vessels; e.g., the middle cerebral artery (MCA), 
with no moving parts. If the sensor is pointed to the general volume location of the 
15 desired blood vessel (e.g., within ± 1 cm.), it will lock to within ± 0.1 mm of the point 
of maximum radial component of blood flow and remain locked in spite of patient 
movement. 

A device of the present invention can remain focused on the selected blood vessels 
regardless of patient movement because it produces and digitally analyzes, in real 
20 time, a 5-dimensional data base composed of signal-return amplitude as a function 
of: 

1. Depth, 

2. Azimuth, 

3. Elevation, 

25 4. Radial component of blood velocity, 
5. Time. 

Since a device of the present invention can automatically locate and lock onto the 
point with the maximum volume of blood having a significant radial velocity, 
unattended continuous blood velocity monitoring is one of its uses. By using the 
30 precise relative location of the point at which lock occurs as a function of depth, a 
device of the present invention can map the network of blood vessels as a 3- 
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dimensional track without the hardware and computational complexity required to 
form a conventional ultrasound image. Using the radial component of velocity along 
with the three-dimensional blood path, a device of the present invention can directly 
compute vector velocity. 

5 A device used in a method of the present invention is a non-mechanical Doppler 
ultrasound-imaging sensor comprising probes, processing electronics, and display. 
Specific choices of probes allow the system to be used for transcranial Doppler 
(TCD), cardiac, dialysis, and other applications. 

10 The present invention may be better understood by reference to the following non- 
limiting Examples, which are provided as exemplary of the invention. The following 
Examples are presented in order to more fully illustrate particular embodiments of the 
invention. They should in no way be construed, however, as limiting the broad scope 
of the invention. 

15 EXAMPLE 1 AN ULTRASOUND DIAGNOSTIC AND MONITORING 

SENSOR WITH REAL-TIME 3-D MAPPING AND TRACKING 
OF BLOOD FLOW 

This embodiment of the present invention has application for medical evaluation and 
20 monitoring multiple locations in the body; however, the transcranial Doppler 
application will be used as an example to describe the invention. 

This invention provides: (1 ) affordable three-dimensional imaging of blood 
flow using a low-profile easily-attached transducer pad, (2) real-time vector 
velocity, and (3) long-term unattended Doppler-ultrasound monitoring in spite 
25 of motion of the patient or pad. None of these three features are possible with 
current ultrasound equipment or technology. 

The pad and associated processor collects and Doppler processes ultrasound 
blood velocity data in a three dimensional region through the use of a planar 
phased array of piezoelectric elements. Through use of unique beamforming 
30 and tracking techniques, the invention locks onto and tracks the points in 
three-dimensional space that produce the locally maximum blood velocity 
signals. The integrated coordinates of points acquired by the accurate 
tracking process is used to form a three-dimensional map of blood vessels 
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and provide a display that can be used to select multiple points of interest for 
expanded data collection and for long term continuous and unattended blood 
flow monitoring. The three dimensional map allows for the calculation of 
vector velocity from measured radial Doppler. 

5 A thinned array (greater than half-wavelength element spacing of the 

transducer array) is used to make the device inexpensive and allow the pad to 
have a low profile (fewer connecting cables for a given spatial resolution). 
The same physical array can also be used to form a broad transmit beam 
encompassing a plurality of narrow receive beams. Initial acquisition of the 

10 blood velocity signal is attained by insonating a large region by defocusing the 
transmit array or by using a small transmitting sub-aperture, for example. The 
computer simultaneously applies numerous sets of delays and/or complex 
weights to the receiver elements in order to form M simultaneous beams. 
With M beams being formed simultaneously, the receiver can dwell M times 

is as long, so as to obtain high S/N and fine Doppler resolution. For an 

embodiment that utilizes a small transmitting sub-aperture, the source of the 
transmitted energy within the array (i.e., the location of the transmitter sub- 
aperture) varies with time in order to lower the temporal average spatial peak 
intensity to prevent skin heating. 

20 The array is thinned without reducing the receiver area by limiting the angular 
field of view. When needed, a map of a larger region is created by moving the 
probe and correlating the sub-images. Once the region of interest has been 
imaged and coarsely mapped, the full transmitter array is focused at a 
particular location on a particular blood vessel for tracking. In the tracking 

25 mode: (1) grating lobes due to array thinning are reduced by using wide 
bandwidth and time delay steering and (2) only three beams are formed at a 
time: sum, azimuth difference, and elevation difference. When the desired 
(i.e. peak) blood velocity signal is not in the output, this is instantly recognized 
(e.g., a monopulse ratio, formed after Doppler filtering, becomes non-zero) 

30 and the array is used to track (slow movement) or re-acquire (fast movement) 
the desired signal. Re-acquisition is achieved by returning to step one to form 
and Doppler-process a plurality of beams in order to select the beam (and the 
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time delay or "range gate") with the most high-Doppler (high blood velocity) 
energy. This is followed by post-Doppler monopulse tracking in azimuth, 
elevation, and range to lock a beam and range gate on to the exact location of 
the peak velocity signal. 

5 In applications such as transcranial Doppler, where angular resolution based 
on wavelength and aperture size is inadequate, fine mapping is achieved, for 
example, by post-Doppler monopulse tracking each range cell of each vessel, 
and recording the coordinates describing the location of the monopulse null. 
With a three-dimensional map available, true vector velocity can be 

10 computed. For accurate vector flow measurement, the monopulse difference 
is computed in a direction orthogonal to the vessel by digitally rotating until a 
line in the azimuth-elevation or C-scan display is parallel to the vessel being 
monitored. 

All current ultrasound devices (including "Doppler color flow mapping" 
15 systems) form images that are limited by their resolution. In some 

applications, such as TCD, the low frequency required for penetration makes 
the azimuth and elevation resolution at the depths of interest larger than the 
vessel diameter. In this invention, as long as (1) a blood vessel or (2) a flow 
region of a given velocity can be resolved by finding a 3-D resolution cell 
20 through which only a single vessel passes, that vessel or flow component can 
then be very accurately located within the cell. Monopulse is merely an 
example of one way to attain such sub-resolution accuracy (SRA). SRA 
allows an extremely accurate map of 3-D flow. 

This invention utilizes post-Doppler, sub-resolution tracking and mapping; it 
25 does Doppler processing first and uses only high Doppler-frequency data. 
This results in extended targets since the active vessels approximate "lines" 
as opposed to "points". In three-dimensional space, these vessels are 
resolved, one from another. At a particular range, the azimuth-elevation axis 
can be rotated so that the "line" becomes a "point" in the azimuth dimension. 
30 That point can then be located by using super-resolution techniques or by 
using a simple technique such as monopulse. 
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Overview of the Embodiment 

The invention is complex because it involves (1) a family of ultrasound 
sensors (for different parts of the body), (2) the interplay of a set of core 
technologies that are unique by themselves, and (3) a number of design 
5 options which represent different ways to implement the invention. To 
facilitate an organizational understanding of this many-faceted invention, we 
precede a description of an overall preferred embodiment with a discussion of 
each of the three topics above. 

The sensors addressed are all two-dimensional (i.e., planar) arrays of 
10 piezoelectric crystals for use in active, non-invasive, instantaneous (or real- 
time), three-dimensional imaging and monitoring of blood flow. While the 
sensors and the techniques for their use apply to all blood vessels in the 
body, the figures and detailed description emphasizes the transcranial 
Doppler (TCD) monitor because that application is most difficult to implement 
15 without all of the components of this invention. The sensors use a unique 

approach to 3-D imaging of blood velocity and blood flow that (1) allows for 

i 

finer image resolution than would otherwise be possible with the same 
hardware complexity (number of input cables and associated electronics) and 
(2) allows for finer accuracy than would ordinarily be possible based on the 
20 resolution. The invention measures and monitors 3-D vector velocity rather 
than merely the radial component of velocity. 

The core technologies that constitute the invention are (1 ) array thinning with 
suppression of ambiguities or grating lobes, (2) post-Doppler sub-resolution 
tracking, and (3) post-Doppler sub-resolution mapping. The invention 

25 encompasses two ways to thin the array (reducing the number of input cables 
and input signals to be processed while maintaining high resolution and 
avoiding ambiguities). The first is bistatic operation; the second is broadband 
operation. In the TCD application, where signal to noise and hence receiver 
array area is of paramount importance, array thinning is possible without 

30 reducing the receiver array area because a relatively small (compared to 
other applications) angular field of view is needed. One particular bistatic 
approach to thinning reduces transmitter area and consequently poses a 
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problem of excessive spatial peak intensity (skin heating) in the TCD 
application. This is solved by a component invention called transmitter 
diversity (which lowers the temporal average of the spatial peak intensity). 
The phase-defocusing bistatic approach and the monostatic or bistatic 
5 broadband approach to thinning all use the entire aperture and hence do not 
require transmitter diversity. 

In the TCD application, the achievable angular resolution is poor, regardless 
of the method of thinning, or whether or not thinning is used. Once a section 
of a blood vessel is resolved from other vessels in Doppler, depth, and two 
10 angles (az and el), Post-Doppler sub-resolution processing locates that 
section to an accuracy that is 10 to 20 times as fine as the resolution. This 
allows for precise tracking and accurate mapping. Tracking provides for the 
possibility of unattended long term monitoring and mapping aids the operator 
in selecting the point or points to be monitored. 

is There are many options available in the design of any member of the family of 
sensors that utilizes any or all of the core technologies that comprise this 
invention, A two-dimensional array is established art that can be designed in 
many ways and can have many sizes and shapes (rectangular, round, etc.). 
Digital beamforming (DBF) is a technique that has been in the engineering 

20 literature (especially radar and sonar) for many years. One medical 
ultrasound DBF patent cites many references, while another describes a 
particular instance of DBF without citing the other patent or any other prior art. 
While planar arrays, DBF, Doppler ultrasound, and color flow imaging are 
prior art, the manner in this specification of using such established 

25 technologies to map, track, measure, and monitor blood flow is unique. 
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The embodiment is a non-invasive, continuous, unattended, volumetric, blood 
vessel tracking, ultrasound monitoring and diagnostic device. It will enable 
unattended and continuous blood velocity measurement and monitoring as 
5 well as 3-dimensional vascular tracking and mapping using an easily 

attached, electronically steered, transducer probe that can be in the form of a 
small pad for monitoring application, when desired. Although the device has 
application to multiple body parts, the cranial application will be used as a 
specific example. The device can, for example: 

10 1 . Measure and continuously monitor blood velocity with a small low-profile 
probe that can be adhered, lightly taped, strapped, banded, or otherwise 
easily attached to the portion of the body where the vascular diagnosis or 
monitoring is required. 

2. Track and maintain focus on up to four desired blood vessels in spite of 
15 movement. 

3. Map 3-D blood flow; e.g., in the Circle of Willis (the central network of 
arteries that feeds the brain). 

4. Perform color velocity imaging and display a 3-D image of blood flow that 
is rotated via track ball or joystick until a desired view is selected. 

20 5. Form and display a choice of projection, slice, or perspective views, 
including (1) a projection on a depth-azimuth plane, a B-scan, or a 
downward-looking perspective, (2) a projection on an azimuth-elevation 
plane, a C-scan, or a forward-looking perspective, or (3) a projection on an 
arbitrary plane, an arbitrary slice, or an arbitrary perspective. 

25 6. Use a track ball and buttons to position circle markers on the points at 
which we wish to measure and monitor vector velocity. 

7. Move the spatial resolution cell being measured along the blood vessel by 
using the track ball to slide the circle marker along the image of the vessel. 

8. Display actual instantaneous and/or average vector velocity and/or 
30 estimated average volume flow. 
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9. Maintain a 3-day history and display average blood velocity versus time 
for each monitored vessel over 14 hours. 

10. Sound an alarm when maximum or minimum velocity is exceeded or when 
emboli count is high. 

5 1 1 .Track, map, and monitor vessels as small as 1 mm in diameter, resolve 
vessels as close as 4 mm apart (for example), and locate them with an 
accuracy of ±0.1 mm. 

The Monitoring Device will allow a person with little training to apply the 
sensor and position it based on an easily understood ultrasound image 

10 display. The unique sensor can continuously monitor artery blood velocity and 
volume flow for early detection of critical events. It will have an extremely low 
profile for easy attachment, and can track selected vessels; e.g., the middle 
cerebral artery (MCA), with no moving parts. If the sensor is pointed to the 
general volume location of the desired artery (e.g., within ± 0.5 cm.), it will lock 

15 to within ± 0.1 mm of the point of maximum radial blood flow and remain 
locked in spite of patient movement. 

The device can remain focused on the selected blood vessels regardless of 
patient movement because it produces and digitally analyzes, in real time, a 
5-dirnensional data base composed of signal-return amplitude as a function 
20 of: 

6. Depth, 2. Azimuth, 3. Elevation, 4. Radial blood velocity, 5. Time. 

Since the device can automatically locate and lock onto the point with the 
maximum volume of blood having a significant radial velocity, unattended 
continuous blood velocity monitoring is one of its uses. By using the precise 
25 relative location of the point at which lock occurs as a function of depth, the 
device can map the network of blood vessels as a 3-dimensional track without 
the hardware and computational complexity required to form a conventional 
ultrasound image. Using radial velocity along with the three-dimensional blood 
path, the device can directly compute vector velocity. 

30 The proposed device is a non-mechanical Doppler ultrasound-imaging sensor 
consisting of probes, processing electronics, and display. Specific choices of 
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probes allow the system to be used for transcranial Doppler (TCD), cardiac, 
dialysis, and other applications. 

Fig. 1 shows the TCD configuration and the initial definition of the display 
screen. The TCD system is comprised of one or two probes attached to the 
5 head with a "telephone operator's band" or a Velcro strap. The interface and 
processing electronics is contained within a small sized computer. A thin 
cable containing 64 micro coax cables attaches the probe to the electronics in 
the computer. When the operator positions the probe on the head the 
Anterior, Middle and Posterior Cerebral Arteries and the Circle of Willis are 

10 imaged on the screen along with other blood vessels. The arteries or vessels 
of interest are selected by viewing the image. The system locks onto the 
blood vessels and tracks their position electronically. A variety of selected 
parameters is presented on the screen; e.g., the velocity, the pulse rate, depth 
of region imaged, gain and power level. Using only one probe the TCD can 

is monitor up to two arteries (vessels) at a time. Presented on the screen are 
dual traces, one for each artery. The blood velocity can be dynamically 
monitored. As shown in Fig. 1 both the current blood velocity (dark traces) 
and any historic trace (lighter color) can be displayed simultaneously. The 
average blood velocity or estimated average flow for each artery is displayed 

20 below the respective velocity trace. The image shows the arteries and the 
channel used for each artery. When two probes are used, the display is split 
showing signals from both of them. Using a different probe (i.e., different size) 
with the same electronics and display, the unit can be used to measure and 
monitor the blood flow in a carotid artery. Similarly, it can be used to perform 

25 this function for dialysis, anesthesia, and in other procedures. 

The sensor is a two dimensional array of transducer elements (piezoelectric 
crystals) that are configured and utilized differently for transmit and receive 
during acquisition. For example, if a square (NxN) array is used, all A/ 2 
elements would receive at the same time, but only a 2x2 sub-aperture would 
30 transmit at any one time. This is illustrated in FIG. 2 for the case of A/=8. The 
array need not be square. Any MxN array may be utilized in this manner. All 
NM received signals (64 in our example) are sampled, digitized, and 
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processed. This can be done, for example, in a desk top or lap top personal 
computer with additional cards for electronics and real-time signal processing 
as illustrated in FIG. 1 and FIG. 3. If the PCI bus in FIG 2 becomes a 
bottleneck for high speed processing, a pipelined or systolic architecture 
5 would be used. Alternatively, the processing can be performed in an 
application specific integrated circuit (ASIC). 

The small (4 element) transmit sub-aperture (FIG. 2 b) produces a broad 
transmit beam that insonates a region containing many receive beams. This 
is schematically illustrated in FIG. 4 for the particular case of a square array 

10 and square elements such as in FIG. 2. Since data is received from each 
element of the array, this data can be combined in a processor (FIG. 3, for 
example) in many different ways to form any number of beams. The 
transmitter is larger than a single array element so that it can provide some 
selectivity and not insonate the grating lobes caused by array thinning 

15 (spacing the array elements more than 14 wavelength apart). The concept is 
illustrated below for a 1 -dimensional array forming a beam that measures only 
one angle. For a two-dimensional array, this represents a horizontal or 
vertical cut through the cluster of beams shown in FIG. 4. FIG. 4 was an 
approximate and conceptual representation of the two-angle (azimuth and 

20 elevation) extension of the single angle case detailed below. 

"Grating lobes" are ambiguities or extra, unwanted, beams caused by using a 
transducer array whose elements are too large and hence too far apart. The 
following analysis illustrates grating lobe suppression for the worst case of 
narrowband signals and phase-shift beam processing. Time delay processing 
25 using wideband signals would be similar, but would further attenuate or 
eliminate grating lobes, resulting in even better performance. 

The next four figures show beam pattern amplitudes plotted against 

x = (d/X)sm0, (1) 

30 



where x represents a normalization for the angle, 0, from which reflected 
acoustic energy arrives. The azimuth (or elevation) angle, 0, is zero in the 
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broadside direction, perpendicular to the transducer array. The width (or 
length) of a transmitter is 2d, where d is the width (or length) of a single 
element of the receiver array. The wavelength of the radiated acoustic wave 
is X = c/f 9 where c is the acoustic propagation velocity (1 540 meters/second in 
5 soft tissue) and /is the acoustic frequency (usually between 1 and 10 
megahertz). FIG. 5a shows the transmitter pattern 

arfx) = sin 2ml 2txx (2) 

10 for the special case of uniform insonation over the 2d-wide transmitter sub- 
aperture being used. 

The receiver pattern is the product of the receiver element pattern and the 
receiver array pattern 

1 5 a R (x) = ajuix) a^x) (3) 

Each of these two component patterns is plotted separately in FIG. 5b. Again 
assuming the special case of a uniform receiver element (and a square 
element in the case of a 2-D array), the receive element pattern is 

20 

aR£(x) ~ sin 7tx / tdc. (4) 

The receiver element pattern is twice as wide as the transmitter pattern 
because the receiver element is half as wide as the transmitter. In the far- 
25 field, i.e., for Ar»L 2 t where r is the range or depth and L is the length of the 
aperture, the receive array pattern steered to the angle 0= 6b is 

flwW-Sw.^^^, (5) 

30 where w n is a weighting to reduce sidelobes and N is the number of elements 
in one dimension. As seen in Figure 5b, equation (5) is periodic in x. The 
peak at x=x 0 (x 0 = 0 in Fig 5) is the desired beam and the others are grating 
lobes. 
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In the near field, when focused at (r 0 ,6b), equation (5) is replaced by the 
slightly better general Fresnel approximation: 



(provided that that the range significantly exceeds the array size, r>L), where 
x = d sin 01 X y as before, and 



Because the receiver aperture is sampled with a spatial period of d t the 
receiver array pattern will be periodic in sin 0, with a period of A/d (equation 5). 
This periodicity means that the array pattern is ambiguous. When the array is 
pointed broadside (0= 0), it will also be pointed at the angle 0= sin' 1 (A/d), for 
example. In terms of the normalized variable, *, the period is unity. Since | sin 
6\ cannot exceed 1 , the variable x is confined to the interval [-d/A, d/A]. 
The conventional element spacing is d = A/2. Thus, in a conventional phased 
array,* is always between -0.5 and +0.5, and hence ambiguities are not 
encountered. In a highly thinned array (d> X) t there will normally be 
ambiguities or grating lobes as illustrated in FIG. 5b. The second grating lobe, 
at x=2 or 0 = sin* 1 (2 A/d), is not real when d does not exceed 2A. 

FIG. 5c shows the two-way pattern. The gating lobe suppression, resulting 
from the choice of a transmitter diameter of D = 2d is valid for all values of d. 
In a two dimensional array, the elements could be rectangular instead of 
square (d x x d y \ and the results would still be valid. Similar results could be 
obtained for an array in which the elements are staggered from row to row 
(and/or column to column). For example, if the receiver array is a "bathroom 
tile" of hexagonal elements, the transmitters could be chosen as sub-arrays 
consisting of an element and its six surrounding neighbors. 

In FIG. 6 the same array is used as in FIG. 5, but the receiver element signals 
are combined with a phase taper that steers the beam to x = 0.2. This is 
approximately (a little less than) the half power point, where a£x) a n (x) = 




(6) 



0/Ar. 



(7) 
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0.707. In FIG. 6c t we see that the grating lobes are not completely 
suppressed, with the largest one at x = -1+.2 = -.8. FIG. 7 shows this in 
decibels. The worst-case grating lobe is attenuated by at least 25 dB, even in 
the stressing case of extremely narrow band operation. A Hanning window 
5 was applied to keep the sidelobes lower than the peak grating lobe. These 
Figures were produced in MATLAB, using the following software (m-file): 

X=-2:l/64:2-l/64; 
p=pi*x+eps; R=sin(p) . /p; 
10 p=2*p; T=sin(p) ./p; 
N=8 

n=0:N-l; 
% XO=0; 

xo=0.2; % is 2 -way 1/2 power 
15 e=exp(j*n' *2*pi* (x-xo) ) ; 
w=hanning(N) ; 
% E= (1/N) *ones(l,N) *e ; 
E=(2/N)*w'*e; 

subplot (311) ; plot (x,abs(T) ) ; 
20 subplot(312) ; plot (x, [abs (R) ; abs (E) ] ) ; 
TRE=abs (T) . *abs (R) . *abs (E) ; 
subplot (313) ; plot(x,TRE); 
figure (2); plot (x, 20*logl0 (TRE) ) ; 
zoom on; 

25 

The dimensions in FIG 4 are representative for a transcranial Doppler 
application of the invention, to provide a specific example. If /= 2 MHz is 
chosen for the center frequency, the wavelength is 0.77 mm. An 8x8 array 
with a width and/or length of 1=1 cm, provides a one dimensional thinning 
30 ratio of 2 d/X = 3.247. For a square array, the total number of elements is 
reduced by a factor of (2 d/X) 2 > 10 from that of a filled array. Even greater 
thinning ratios are possible. Even if d/X is kept less than 2 to avoid a second 
grating lobe (at x=2), complexity reductions up to a factor of 16 are possible. 
For the 1 cm array at 2 MHz, the hyperfocal distance (where the 3 dB focal 
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region extends to infinity) is L 2 /4A = 3.25 cm. Thus, a fixed focus probe 
suffices for this application. However, since the simultaneous formation of 
multiple receive beams is conveniently performed digitally, dynamic focus on 
receive is easily accomplished. The quadratic phase distribution across the 
5 elements required to focus in depth is simply added to the linear phase 
distributions required to steer the beams. 

FIG. 8a shows the product of the transmitter pattern (FIG. 5a or 6a) and the 
receiver element pattern. FIG. 8b plots the element patterns for a set of five 
10 beams steered to x = -.2, -.1 , 0, ,1, and .2. This set of five receive beams 
shows grating lobes of the thinned array. FIG. 8c shows the set of resulting 2- 
way patterns obtained by multiplying the patterns in FIG. 8b by the function 
plotted in FIG. 8a. Here, the grating lobes are suppressed. This represents a 
horizontal or vertical cut through the cluster of beams in FIG. 4. 

15 

Using the configuration described above, the cluster of beams in Figs. 4 and 
8c is used to approximately locate the desired point for collecting the blood 
velocity signal. For example the output of each beam in the cluster would be 
Doppler processed by performing an FFT or equivalent transformation on a 

20 sequence of pulse returns. The pulse repetition frequency (PRF) would 
typically be less than or equal to 9 kHz to unambiguously achieve a depth of 
8.5 cm for the TCD application. In order to obtain a velocity resolution as fine 
as Av = 1 cm per second (to distinguish brain death), a dwell of duration 7 = 
JU(2 Av) = 38.5 ms, or 347 pulses at 9kHz, is desired. For efficient FFT 

25 processing, the number of pulses used would be zero filled to a power of 2 
such as 512. 

The example shown in Figs. 2 through 8 was an 8 by 8 receiver array forming 
a 5 by 5 cluster of beams. This is an example of an approximate rule of 
30 thumb for this invention, that an N element linear array is recommended for 
use in producing N/2 +1 beams for N even and [7/+1]/2 beams for //odd. 
Thus, a 16 by 10 element rectangular array would preferably be used to form 
a 9 by 6 cluster of beams, though the actual number of beams formed is 
arbitrary. This recommended number of beams is derived below. 
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If an N elements were used to form orthogonal beams, e.g., by an TV-point 
FFT, then there would be N beams in a 180° angular region, from -90° to 
+90°, corresponding to -1 <w<1 , where u = sin 0. In conventional phased array 

5 ultrasound, a 128 (=N) element array is used to produce 256 (=2/V) lines 
(sequentially scanned beams) in a 90° angular region from -45° to 45°, 
corresponding to -.707<w<.707. If the array is filled, then x = ul2 (Equation 1) 
and 2N beams are conventionally formed in \x\<^2 14. When we thin the array, 
we prefer to have |x|<0.2 = 1/5 (the 3 dB point of the curve in Fig. 7a). The 

io number of beams in that region, for the same beam density as used in current 
practice, is given by 

Recommended No. of beams = (M5)N -f (V2 1 A) = 2 V2 N/5 * 0.5657 N. 

is The beams are formed digitally, using software on a personal computer or 
using digital signal processing hardware to implement equations such as 
Equation 5 or 6. The electronic interface between the probe and the 
processor is diagrammed in Fig. 9. This figure illustrates the case of signals 
from 64 elements being connected to a single A/D converter, and power being 

20 applied to sets of four elements. The use of a separate A/D converter for 
every received channel, for example, is another possible implementation of 
this invention. 

A conventional, half-wavelength spaced, monostatic, phased array could 
25 sequentially search a region of interest, but it would require far more elements 
and would thus be far more costly. Using the array differently in transmit and 
receive, not only allows for the formation of multiple beams; it also enables 
the use of the angular pattern of the transmitter to suppress receiver grating 
lobes. This allows for a "thinned" array (elements spaced less than a half 
30 wavelength apart). Because receive beams are formed only in a limited 

angular region, a wide-angle receiver element pattern (which usually implies a 
small element) is not required. In fact, the size of the receiver element can be 
as large as the element spacing. Thus the receiver array is "thinned" only in 
the sense that the element spacing exceeds a half wavelength. Since the 
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element size also exceeds a half wavelength, the array area is not reduced. 
It is thinned only in terms of number of elements, not in terms of receiver area. 
Consequently, there is no reduction in signal-to-noise ratio, nor a requirement 
for increased transmitter power. 

5 

A monostatic array would transmit from the full aperture, scanning the 
transmitted beam over the region being examined. The "bistatic" array of this 
invention transmits from a sub-aperture to insonate multiple receive beam 
positions simultaneously. Since there is an FDA limit to spatial peak, 

10 temporal average, intensity (l spta ), there may be a danger of exceeding this 
limit at the transducer surface, creating a danger of burning the skin. This 
potential danger is eliminated by using a different transmit sub-aperture for 
each coherent dwell or burst of pulses. This transmitter diversity technique 
spreads the temporal average intensity over the face of the array, reducing 

15 I spta to what it would be if the entire array were used at once. 

For the particular implementation pictured in FIG. 9, an A/D converter is 
multiplexed amongst the 64 elements. The signal spectrum at any of these 
elements is centered at/ 0 = 2 MHz, as shown in FIG. 10a. This is a real signal 

20 with a spectrum that is symmetric about/= 0. This analog signal is bandpass 
filtered (BPF) to insure that there is little power outside of a 444 kHz band 
centered at 2 MHz. If a 51 2/9 =56.889 MHz A/D converter is used, each 
receive channel is sampled at/, = 888.9 kHz, giving rise to a real sampled 
signal with a spectrum as shown in FIG. 10b. A processing element such as 

25 a field programmable gate array (FPGA) is used to shift the frequency by/,/4 
(FIG. 10c) by "multiplying" by quarter cycle samples of sinusoids (which are 
zeros and ones). The same FPGA also digitally filters (or Hilbert transforms) 
the complex signal to decimate its sampling rate by a factor of two. The 
spectrum of the decimated digital low-pass signal is shown in FIG, 10 d. 

30 The signal sent to the processor from each element has the spectrum shown 
in Fig. 9d, and consists of r=f s /2 complex samples per second. The total data 
rate into the processor is approximately 57 megabytes per second. For non- 
real-time operation, tens of seconds of data at a time will be collected in 
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system memory and then transferred to hard disk. For real-time monopulse 
tracking, only three beams are formed, so that the data rate is reduced to 3 x 
0.8889 = 2.67Mbytes, or 5.33 Mbytes allowing for bit growth. 

5 The transmitted pulses are sent to a group of four elements. The particular 
embodiment shown in Figure 9 uses diodes to block the received signals and 
prevent mutual coupling between the four receive elements. After a coherent 
pulse train (or pulse burst used for Doppler processing), the waveform is 
switched to another set of 4 elements for the next burst. A separate power 

10 amplifier is associated with each of the 1 6 sets of elements so that the 
switching can be accomplished at low power. 

One embodiment of sub-resolution tracking (i.e., tracking and locating blood 
flow to a small fraction of a spatial resolution cell) is "Monopulse". Monopulse 
tracking is performed as follows. A particular set of complex weights are 

is applied to the set of received signals (64 in the example of FIG 2) to steer a 
beam at the middle cerebral artery, for example. The phase taper across the 
array defines the steering direction and the amplitude taper (called a window 
in radar and a shading in sonar) is used to provide low sidelobes for high 
dynamic range. The beam output (a linear combination of the signals) is 

20 range gated (time delay corresponding to the desired depth) and the range- 
gated / beam-formed output from a sequence of transmitted pulses is then 
Fourier transformed to obtain a plot of amplitude versus Doppler frequency. 
The receive beam is steered digitally to the point that produces the maximum 
amplitude at high Doppler frequencies. 

25 Since the measured data at each element is stored, the digital processor can 
apply more than one set of weights at a time, forming more than one beam. 
For software monopulse the processor will form three beams, all in the same 
direction. All three beams may have the same phases applied to the element 
signals; but the amplitudes will differ. The beam called Sum has all positive 

30 amplitudes, with the larger weights applied to the central elements. This 
forms a fairly broad beam. The b am called Az for "azimuth difference beam" 
has large positive weights on the rightmost elements and large negative 
weights on the leftmost elements (or vice versa). The beam called £/for 
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elevation difference has large positive weights on the top-most elements and 
large negative weights on the bottom-most elements. A correctly pointed 
beam would have Az = El- 0, and Sum would be maximized. 

The ratio of the peak Doppler amplitude outputs: Az/Sum, is a precise measure 
of the azimuth pointing error and the corresponding ratio El/Sum measures the 
elevation pointing error. The digital steering phase taper is thus corrected 
with data from a single burst of pulses. The duration of the pulse burst is the 
reciprocal of the medically required Doppler resolution (usually corresponding 
to the minimum blood velocity that can support life). Without techniques such 
as those described in this specification, a sequence of at least four additional 
Doppler dwells or pulse bursts would be required (above, below, to the right, 
and to the left) in a hunt and seek method to find the correct (maximum peak 
Doppler Amplitude) beam. With monopulse, the correction is very precise (to 
within ±0.1 mm of the point of maximum peak Doppler amplitude) and virtually 
instantaneous. For the bistatic digitally beamformed sensor, the original data 
exists in computer memory. Hence, whenever the Doppler processed 
monopulse differences are non zero, the same data set could even be re- 
processed to form a correctly pointed beam. A slower processor would 
merely process the next burst correctly. 

A "front view" perspective display or a C scan display (azimuth horizontal and 
elevation vertical) of the blood flow map at the desired range will allow 
someone to aim the transducer probe or pad at the desired point (highest 
amplitude for high Doppler), so that the desired point is initially within the 
center beam. The receiver array is then steered electronically so that the 
monopulse differences are zero and hence the central beam is precisely 
aimed at the desired point. Slight motions are corrected using monopulse and 
large motions are corrected by again forming all beams to re-acquire the peak 
signal. All corrections are made entirely electronically, in the data processing 
or digital beamforming. A narrow receiver beam will always be precisely 
pointed (to within a tenth or 20 th of the receiver beamwidth) as long as the 
desired point remains within the much larger region covered by the transmitter 
(FIG. 4). 
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True vector velocity is computed from the blood vessel map and the radial 
velocity measured from the pulse Doppler dwell. A map, far more accurate 
than that attainable with the available angular resolution is attained as follows. 
The low-resolution map is used to locate a vessel of interest and a beam is 

5 locked on it at a fixed range, using azimuth and elevation monopulse. The 
coordinates of the point at which lock occurs is recorded. The range is then 
changed slightly, another lock (on the same vessel) is obtained, and the 
coordinates are recorded. In this manner, the vessel is mapped far more 
accurately than would be predicted from the available image resolution. All 

10 vessels within the field of view of the probe are similarly mapped. By moving 
the probe angle slightly, another region can be mapped in the same manner. 
Several such maps can be correlated over the region of pair-wise overlap and 
converted to a common coordinate system. In this manner a larger region is 
mapped and displayed than that of the current field of view. The current field 

15 of view would be highlighted, outlined, or presented as a color flow map. 
Points to be monitored in the current field are then selected by moving a 
cursor along the display (point and click). The selected points are Doppler 
processed and tracked using three-dimensional monopulse. While Doppler 
measurements provide only the radial component of velocity, the accurate 

20 blood vessel map provides the exact three-dimensional orientation of the 
vessel at the point being monitored. The measured radial velocity is divided 
by the projection of a unit vector representing the vessel at the monitored 
point onto the transducer line of sight. This gives the magnitude of the true 
vector blood velocity. 

25 Sub-resolution mapping accuracy is attainable if (1 ) the range-azimuth- 
elevation-Doppler resolution cell being examined encompasses only a single 
blood vessel, and (2) "azimuth" monopulse is performed with the usually 
vertical e-axis tilted so that the orientation of the vessel in the spatial 
resolution cell being processed is parallel to the e-r plane ("azimuth" is 

30 constant). 

The user will ascertain from the display, that the resolution cell being 
monitored contains only a single vessel, and would rotate the 3-D blood- 
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vessel map to a C-scan aspect (elevation up and azimuth to the right). A 
vertical mark will appear in the display, within the resolution circle, to signify 
the orientation of the monopulse axis. This axis (parallel to the line separating 
the positively and negatively weighted array elements) can then be oriented 
5 so that the mark is aligned with the blood vessel in either of two ways. The 
probe can be physically twisted (rotated about the line of sight), or it can be 
electronically rotated via digital processing because the weights are applied 
digitally. 

FIG. 1 1 illustrates the segment of a vessel in a single resolution cell, after 
10 rotation. The resolution cell shown is not a cube because the range resolution 
will typically be finer than the cross-range resolution. The illustrated circular 
cylinder represents blood cells in a vessel reflecting energy at a fixed Doppler 
frequency. These represent a cylindrical annulus of blood cells, at a constant 
distance from the vessel wall, moving with approximately the same velocity. 
15 In the single resolution cell of Fig. 1 1 , the return at the highest Doppler would 
represent a line in three-dimensional space (the axis of the vessel) and hence 
a point on the azimuth axis after rotation. When applied to the highest 
Doppler output, the Sum beam would have broad peak at zero azimuth (a = 0) 
and the monopulse ratio, r-Az/Sum, will be a linear function of the azimuth 
20 angle to which the array is phase steered: 

r(a)=k a. 

This result can be attained by applying the same phase across the aperture 
for the Az and Sum beams, but using the derivative of the Sum beam amplitude 
weights with respect to x and y respectively for the Az and El aperture weights. 

25 Other Embodiments 

If the wide transmit beam (for search and acquisition) is created by using a 
quadratic phase curvature instead the scheme of Fig. 2b, transmitter diversity 
may not be needed. Furthermore the manner of controlling grating lobes in 
Fig.1 and Figs.5 - 8 is only one of many. Using a wider bandwidth and time- 
30 delay steering can also reduce grating lobes. 
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EXAMPLE II ULTRASOUND MEASUREMENT OF BLOOD VOLUME 

FLOW 

As explained above, current ultrasound Doppler devices measure radial velocity. 
Several methods now exist for 3-D ultrasound imaging, such as those involving 

5 transducer motion. A three-dimensional image with Doppler allows for the 

measurement of vector velocity. Example I above provides measurement and long 
term monitoring of three-dimensional vector velocity. If the resolution of a color flow 
Doppler image is sufficient to provide an estimate of the inside diameter of the blood 
vessel, then measurement of volume blood flow becomes practical. Present ly 

10 available ultrasound imaging devices have either low resolution or they only produce 
a two-dimensional image. The present invention combines vector velocity 
information (such as attained as explained in Example I above) with additional 
information to obtain volume flow. The additional information is the inside diameter 
of the vessel under examination, the blood velocity profile across the vessel, or the 

15 vector velocity as a function of time and position (i.e., the velocity field). This 

additional information can be obtained from a high-resolution radial-Doppler or color 
flow image or from external data such as a high-resolution MRI image. 

A two-dimensional array of piezoelectric elements, or some other means, is used to 
image blood flow in a three dimensional region. A particular point on a particular 
20 vessel is selected and the vector representing the orientation of the vessel is noted. 
The radial velocity divided by the cosine of the angle made by the vessel with the line 
of sight at the measurement point is the magnitude of the vector velocity. That 
number integrated over the vessel cross section would give the volume flow in 
volume per unit time or milliliters per minute, for example. 

25 Figure 12 shows a circular cylinder representing blood cells in a vessel moving at a 
particular velocity and thus reflecting energy at a specific Doppler frequency. The 
figure assumes that methods such as those in the referenced invention, for example, 
have been used to measure the 3-D orientation of the vessel so that the vector 
velocity can be calculated and the azimuth axis can be defined to be perpendicular to 

30 the vessel. 

The simplest way to estimate volume flow is to measure the vessel diameter, rf, (or 
radius d/2), calculate the cross-sectional area, A = n(d/2) 2 , and multiply by the 
average velocity. A more accurate way is to integrate the velocity as a function of 
position, over the cross-section. The velocity is a function of the radius, a, of the 
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10 



15 



20 



cylinder depicted in Fig. 12. If a is the distance from the cylinder to its axis, and v 
(a) is the velocity function, then the volume flow is 



Equation (1) assumes a circular cross-section of constant radius, p=d/2. It is a 
special case of the more general polar coordinate integration: 



The velocity function is determined by determining the diameter (and hence the 
radius) of the cylinder corresponding to each velocity. 

For example, a 1.5-cm diameter Doppler ultrasound transducer array operating at 10 
MHz will be oriented with the length or azimuth direction perpendicular to the vessel 
to produce a B-scan (depth-azimuth) image. At a depth of approximately 10 mm, the 
cross range resolution is 0.1 mm. If the vessel diameter is 1 mm, the diameter can 
be measured with an accuracy of ± 5%. The area of the vessel is thus known to an 
accuracy of 10%. Since the average vector velocity can be measured extremely 
accurately, the volume flow is also accurate to ± 10%. The best accuracy is attained 
by measuring the azimuth extent corresponding to various velocities and then 
numerically evaluating equation (7) or (8). Naturally, a skilled artisan can readily 
program a processor to solve these equations, and calculate blood flow volume using 
routine programming techniques. 

Since the autocorrelation function (pulse-to-pulse, at a fixed range) and the Doppler 
Power Spectrum form a Fourier pair, the total power can be obtained either as the 
autocorrelation function at zero lag or the integral of the Doppler Power Spectrum 
(Spectral Density) over all Doppler frequencies. Since radial velocity is proportional 
to Doppler frequency, the mean velocity can be obtained from the autocorrelation 
function as shown below: 




(7) 



0 




(8) 
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J- lS(G>)e J °»dco=)s d (f)e J2 * r df, 



2* t 

— oo 

hence 

oo 

= P d = total Doppler power 

-CO 

and 

00 

R» (0 = & a (r) = jt/2*S, (/)>> 2 *4f , 

-oo 

leading to 

00 

K^) = j\2nfS d {f)df. 

-oo 

Thus 

-00 

where the Doppler frequency and its mean (expected value) are related to the 
radial blood velocity and its mean by 

c 

Hence 



which is used in the autocorrelation method of color-flow blood-velocity imaging. We 
note that if we do not normalize by dividing by the total Doppler power, we obtain a 
power-velocity product that indicates the volume flow rate. This is due to the fact that 
power is directly proportional to area [see Reference 1]. 

Since all velocity vectors are parallel at the narrowest point (the vena contracta), flow 
at that particular point can be considered as non-turbulent, even though severe 
turbulence exists before and after. Reference 1 shows that regurgitant blood flow 
through the mitral heart valve can be quantitatively measured by observing the 
Doppler spectrum at that point and using the power-velocity-integral relation below. 
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In terms of the velocity power spectrum, P{v) = (2/b/c) S d (0, Reference 1 shows that 
the blood vessel area in a "slice" perpendicular to the line of sight is directly 
proportional to the total Doppler power (the total power at the output of the high-pass 
5 wall filter). 

where A 0 and P 0 are the known area and measured power in a narrow beam, smaller 
10 than the vessel. If the blood flow makes an angle 9 with the line of sight, the area, 
and hence power, is increased by the factor 1/cos6. This offsets the fact that only the 
radial component of velocity is measured, so that the power velocity integral provides 
true volume flow: 

15 Q = dQ/dt = %jvP(v)dv 

In Reference 1, P was measured with the same probe as P 0 by masking the outside 
of the aperture in order to create a wider beam. With our 2-D phased array, we 
would merely turn off or ignore some of the outer elements. More importantly, we 
20 can use the 3-D image to precisely locate the vena contracta, and lock on to it using 
monopulse. We can even monitor the valve during a stress test, while the patient is 
on a treadmill. 

We note here, that there are several ways to measure the volume flow rate. 

25 Reference 1 uses the fact that it is proportional to the integral of the product of the 
velocity and the power per unit velocity, as in the last equation. An other way is to 
recognize that it is equal to the product of the average radial velocity and the total 
projected area that is, in turn, proportional to total Doppler power. Since the total 
Doppler power is used in the denominator of the autocorrelation-based color-flow 

30 velocity map, volume flow rate can be obtained by merely not dividing by the total 
power. If the P pulse return (after MTI or Doppler high-pass or wall filtering) is 



35 



the volume flow rate is proportional to 



f=l ,2, 
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N 
i=2 

The normalization (denominator) that is used to convert this last quantity to mean 
velocity can be 



I 

that is based on a derivation in Reference 2, or a simple power estimate, such as 

2>, 2 +>? 

I 

10 

The point we wish to make here is that by not dividing by a power estimate to obtain 
radial velocity, we obtain volume flow. Current ultrasound Doppler imaging systems 
compute the mean velocity as a ratio, 

15 E(v) = F/P d , 

and display it as a color flow image. Newer imaging systems [2] also display total 
Doppler power (at the output of the wall filter), P d . By not dividing the color flow 
image by P d , we can also display the true volume flow, dQ/dt This is because the 
20 numerator, 

F = P d *E(v), 

is the power-velocity-integral that is directly proportional to the volume flow. 

25 

Determination of the scale factor, AJP^A/P^ dQldt/F, that must multiply F to 
obtain volume flow requires further comment 

Ao is the area of a reference beam. In [1], A 0 is smaller than the blood flow area. We 
30 will describe three normalization approaches. 

1 . Use a single transmit beam, wider than the vessel, and two simultaneous 
receive beams. One receive beam (the measurement beam) is the same as 
the transmit beam and the other (the reference beam) is smaller than the 
35 vessel. 
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2. Use two (sequential or multiplexed) two-way (transmit and receive) beams. 
One (the measurement beam) is wider than the vessel and the other (the 
reference beam is smaller than the vessel. 
3. Use two (sequential or multiplexed) two-way (transmit and receive) beams. 
5 Both are wider than the vessel and the measurement beam is wider than the 

reference beam. 

Let A 0 be the known area of the reference beam, let P 0 and P 1 be the measured 
received power in the reference and measurement beams. In case 1, the transmit 
10 power density is the same for measurement and reference. The receive power is 
proportional to area. If the area of the vessel (in a slice perpendicular to the line of 
sight) is A, it follows that 

A/Aq = Pt/Po. 

15 In cases 2 and 3 f the transmit power density is greater in the reference beam than in 
the measurement beam, but by a known factor. In all three cases, the power 
received in the measurement beam is proportional to the vessel area. In case 2, the 
received reference power also varies with vessel size, but at a different rate than in 
the measurement beam. With proper calibration, correct measurements can be 

20 attained in all three cases. 

[1]. T. Buck, Et al, "Flow Quantification in Valvular Heart Disease Based on the 
Integral of Backscattered Acoustic Power Using Doppler Ultrasound," Proc. IEEE, 
vol.88, no.3, pp.307-330, March 2000. 

25 

[2]. K. Ferrara and G DeAngelis, "Color Flow Mapping", Ultrasound in Medicine and 
Biology, vol.23, no.2, pp.321 -345, March 1997. 

EXAMPLE III 3-D DOPPLER ULTRASOUND BLOOD FLOW MONITOR 

30 WITH ENHANCED FIELD AND SENSITIVITY 

This example sets forth an ultrasound Doppler device and method that enables non- 
invasive diagnosis (the conventional role of ultrasound systems), and also non- 
invasive unattended and continuous monitoring of vascular blood flow for medical 
applications. In particular, the embodiment of the present invention set forth in this 
35 example provides: (1) affordable three-dimensional imaging of blood flow using a 
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low-profile easily-attached transducer pad, (2) real-time vector velocity, and (3) long- 
term unattended Doppler-ultrasound monitoring in spite of motion of the patient or 
pad. None of these three features are possible with current ultrasound equipment or 
technology. 

5 The pad and associated processor collects and Doppler processes ultrasound blood 
velocity data in a three-dimensional region through the use of a two-dimensional 
phased array of piezoelectric elements on a planar, cylindrical, or spherical surface. 
Through use of unique beamforming and tracking techniques described herein, the 
present invention locks onto and tracks the points in three-dimensional space that 

10 produce the locally maximum blood velocity signals. The integrated coordinates of 
points acquired by the accurate tracking process is used to form a three-dimensional 
map of blood vessels and provide a display that can be used to select multiple points 
of interest for expanded data collection and for long term continuous and unattended 
blood flow monitoring. The three dimensional map allows for the calculation of vector 

15 velocity from measured radial Doppler. 

A thinned array (greater than half-wavelength element spacing of the transducer 
array) is used to make a device of the present invention inexpensive and allow the 
pad to have a low profile (fewer connecting cables for a given spatial resolution). 
The array is thinned without reducing the receiver area by limiting the angular field of 

20 view. Grating lobes due to array thinning can be reduced by using wide bandwidth 
and time delay steering. The array, or portions of the array, is used to sequentially 
insonate the beam positions. Once the region of interest has been imaged and 
coarsely mapped, the array is focused at a particular location on a particular blood 
vessel for measurement and tracking. Selection of the point or points to be 

25 measured and tracked can be based on information obtained via mapping and may 
be user guided or fully automatic. Selection can be based, for example, on peak 
response within a range of Doppler frequencies at or near an approximate location. 

In the tracking mode a few receiver beams are formed at a time: sum, azimuth 
difference, elevation difference, and perhaps, additional difference beams, at angles 
30 other than azimuth ( = 0 degrees) and elevation (= 90 degrees). Monopulse is 
applied at angles other than 0 and 90 degrees (for example 0, 45, 90, and 1 35 
degrees) in order to locate a vessel in a direction perpendicular to the vessel. When 
the desired (i.e. peak) blood velocity signal is not in the output, this is instantly 
recognized (e.g., a monopulse ratio, formed after Doppler filtering, becomes non- 
35 zero) and the array is used to track (slow movement) or re-acquire (fast movement) 
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the desired signal. Re-acquisition is achieved by returning to step one to form and 
Doppler-process a plurality of beams in order to select the beam (and the time delay 
or "range gate") with the most high-DoppIer (high blood velocity) energy. This is 
followed by post-Doppler monopulse tracking to lock a beam and range gate on to 
5 the exact location of the peak velocity signal. In applications such as transcranial 
Doppler, where angular resolution based on wavelength and aperture size is 
inadequate, fine mapping is achieved, for example, by post-Doppler monopulse 
tracking each range cell of each vessel, and recording the coordinates and 
monopulse-pair angle describing the location and orientation of the monopulse null. 

10 With a three-dimensional map available, true vector velocity can be computed. For 
accurate vector flow measurement, the monopulse difference is computed in a 
direction orthogonal to the vessel by digitally rotating until a line in the azimuth- 
elevation or C-scan display is parallel to the vessel being monitored. The aperture is 
more easily rotated in software (as opposed to physically rotating the transducer 

15 array) if the aperture is approximately circular (or eliptical) rather than square (or 
rectangular). Also, lower sidelobes result by removing elements from the four 
comers of a square or rectangular array in order to make the array an octagon. 

All currently available ultrasound devices (including "Doppler color flow mapping" 
systems) form images that are limited by their resolution. In some applications, such 

20 as TCD, the low frequency required for penetration of the skull makes the azimuth 
and elevation resolution at the depths of interest larger than the vessel diameter. In 
this invention, as long as (1) a blood vessel or (2) a flow region of a given velocity 
can be resolved by finding a 3-D resolution cell through which only a single vessel 
passes, that vessel or flow component can then be very accurately located within the 

25 cell. Monopulse is merely an example of one way to attain such sub-resolution 
accuracy (SRA). Other methods involve "super-resolution" or "parametric" 
techniques used in "modern spectral estimation", including the MUSIC algorithm and 
autoregressive modeling, for example. SRA allows an extremely accurate map of 3- 
Dflow. 

30 This invention utilizes post-Doppler, sub-resolution tracking and mapping; it does 
Doppler processing first and uses only high Doppler-frequency data. This results in 
extended targets since the active vessels approximate "lines" as opposed to "points". 
In three-dimensional space, these vessels are resolved, one from another. At a 
particular range, the monopulse angle axis can be rotated (in the azimuth-elevation 

35 plane) so that the "line" becomes a "point" in the monopulse angle direction. That 



WO 00/72756 



PCT/US00/14691 



55 

point can then be located by using super-resolution techniques or by using a simple 
technique such as monopulse. By making many such measurements an accurate 3- 
D map of the blood vessels results. 

Methods for extending the angular field of view of the thinned array (that is limited by 
5 grating lobes) include (1) using multiple panels of transducers with multiplexed 
processing channels, (2) convex V-shaped transducer panels, (3) cylindrical shaped 
transducer panel, (4) spherical shaped transducer panel, or (5) negative ultrasound 
lens. If needed, moving the probe and correlating the sub-images can create a map 
of an even larger region. 

10 Active digital beamforming can be utilized, but the implementation depends on a 
choice to be made between wideband and narrowband implementations. If 
emphasis is on high resolution mapping of the blood vessels, then a wide bandwidth 
(e.g., 50% of the nominal frequency) is used for fine range resolution. If emphasis is 
on Doppler spectral analysis, measurement, and monitoring, the map is only a tool. 

15 In this case, a narrowband, low cost, low range-resolution, high sensitivity 

implementation might be preferred. A wideband implementation would benefit in 
performance (higher resolution, wider field of view, and reduced grating lobes) using 
time-delay steering while a narrowband implementation would benefit in cost using 
phase-shift steering. The invention can thus be described in terms of two preferred 

20 implementations. 

In a wideband implementation, time delay steering can be implemented digitally for 
both transmit and receive by over-sampling and digitally delaying in discrete sample 
intervals. In a narrowband implementation, (1) phase steering can be implemented 
digitally (digital beamforming) for both transmit and receive, and (2) bandpass 
25 sampling (sampling at a rate lower than the signal frequency) can be employed with 
digital down-conversion and filtering. 

Overview of this Embodiment. 

This embodiment of the present invention involves (1) a family of ultrasound sensors, 
(2) the interplay of a set of core technologies that are unique by themselves, and (3) 
30 a number of design options which represent different ways to implement the 
invention. To facilitate an organizational understanding of this many-faceted 
invention, a discussion of each of the three topics above follows. 

The sensors addressed are all two-dimensional (i.e., planar or on the surface of a 
convex shape such as a section of a cylinder) arrays of piezoelectric crystals for use 



WO O0A72756 



PCTYUS00/14691 



56 

in active, non-invasive, instantaneous (or real-time), three-dimensional imaging and 
monitoring of blood flow. While the sensors and the techniques for their use apply to 
all blood vessels in the body, the figures and detailed description emphasizes the 
transcranial Doppler (TCD) monitor method as a nonlimiting example. The method 
5 of the present invention utilizes a new, useful and unobvious approach to 3-D 

imaging of blood velocity and blood flow that (1) allows for finer image resolution than 
would otherwise be possible with the same hardware complexity (number of input 
cables and associated electronics) and (2) allows for finer accuracy than would 
ordinarily be possible based on the resolution. The invention measures and monitors 
10 3-D vector velocity rather than merely the radial component of velocity. 

The core technologies that constitute the invention are (1) array thinning with large 
elements and limited scanning, (2) array shapes to reduce peak sidelobes and 
extend the field of coverage, (3) post-Doppler sub-resolution tracking, (4) post- 
Doppler sub-resolution mapping, (5) additional methods for maximizing the angular 

15 field of view, and (6) various digital beamforming procedures for implementing the 
mapping, tracking, and measurement processes. The invention encompasses array 
thinning, where the separation between array elements is significantly larger than half 
the wavelength. This reduces the number of input cables and input signals to be 
processed while maintaining high resolution and sensitivity and avoiding ambiguities. 

20 In the TCD application, where signal to noise and hence receiver array area is of 
paramount importance, array thinning is possible without reducing the receiver array 
area because a relatively small (compared to other applications) angular field of view 
is needed. 

Thinning with full aperture area imposes limitations on the angular field of view. 

25 Methods for expanding the field of view include using more elements than are active 
at any one time. For example, if the electronics is switched between two identical 
panels, the cross-range field of view at any depth is increased by the size of the 
panel. If the panels are pointed in slightly different directions so that overlapping or 
redundant beams are avoided, the field of view is doubled. A generalization of this 

30 approach involves the use of an array on a cylindrical or spherical surface. 

In the TCD application, the achievable angular resolution is poor, regardless of the 
method of thinning, or whether or not thinning is used. Once a section of a blood 
vessel is resolved from other vessels in Doppler, depth, and two angles (az and el), 
Post-Doppler sub-resolution processing locates that section to an accuracy that is 
35 one-tenth to one-twentieth of the resolution. This allows for precis tracking and 
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accurate mapping. Tracking provides for the possibility of unattended long term 
monitoring and mapping aids the operator in selecting the point or points to be 
monitored. 

One of ordinary skill in the art will readily recognize that there are many options 
5 available in the design of any member of the family of sensors that utilizes any or all 
of the core technologies that comprise this invention, all of which are encompassed 
by the present invention. A two-dimensional array is established art that can be 
designed in many ways and can have many sizes and shapes (rectangular, round, 
etc.). 

10 As with other nonlimiting embodiments of the present invention set forth above, this 
embodiment is a non-invasive, continuous, unattended, volumetric, blood vessel 
tracking, ultrasound monitoring and diagnostic device for blood flow. It will enable 
unattended and continuous blood velocity measurement and monitoring as well as 3- 
dimensional vascular tracking and mapping using an easily attached, electronically 

15 steered, transducer probe that can be in the form of a small pad for monitoring 
application, when desired. Although a device of the present invention has 
applications with blood vessels in any part of the body, the cranial application will be 
used as a specific example. A device of the present invention can, for example: 

1. Measure and continuously monitor blood velocity with a small low-profile 
20 probe that can be adhered, lightly taped, strapped, banded, or otherwise easily 

attached to the portion of the body where the vascular diagnosis or monitoring is 
required. 

2. Track and maintain focus on multiple desired blood vessels in spite of 
movement. 

25 3. Map 3-D blood flow; e.g., in the Circle of Willis (the central network of arteries 
that feeds the brain) or other critical vessels in the cranial volume. 

4. Perform color velocity imaging and display a 3-D image of blood flow that is 
rotated via track ball or joystick until a desired view is selected. 

5. Form and display a choice of projection, slice, or perspective views, including 
30 (1) a projection on a depth-azimuth plane, a B-scan, or a downward-looking 

perspective, (2) a projection on an azimuth-elevation plan , a C-scan, or a forward- 
looking perspective, or (3) a projection on an arbitrary plane, an arbitrary slice, or an 
arbitrary perspective. 
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6. Use a track ball and buttons to position circle markers on the points were 
measurement or monitoring of vector velocity is desired. 

7. Move the track location along the blood vessel by using the track ball to slide 
the circle marker along the image of the vessel. 

5 8. Display actual instantaneous and/or average vector velocity, estimated 
average volume flow, and/or Doppler spectral distribution. 

9. Maintain a multi-day history and display average blood velocity versus time 
for each monitored vessel over many hours. 

10. Sound an alarm when maximum or minimum velocity is exceeded or when 
10 emboli count is high; and maintain a log of emboli detected. 

1 1 . Track, map, and monitor small vessels (e.g., 1 mm in diameter), resolve 
vessels as close as 4 mm apart (for example), and locate them with an accuracy of 
±0.1 mm, for example. 

This embodiment of the present invention will allow a person with little training to 
15 apply the sensor and position it based on an easily understood ultrasound image 
display. The unique sensor can continuously monitor artery blood velocity and 
volume flow for early detection of critical events. It will have an extremely low profile 
for easy attachment, and can track selected vessels; e.g., the middle cerebral artery 
(MCA), with no moving parts. If the sensor is pointed to the general volume location 
20 of the desired blood vessel (e.g., within ± 1 cm.), it will lock to within ± 0.1 mm of the 
point of maximum radial component of blood flow and remain locked in spite of 
patient movement 

A device of the present invention can remain focused on the selected blood vessels 
regardless of patient movement because it produces and digitally analyzes, in real 
25 time, a 5-dimensional data base composed of signal-return amplitude as a function 
of: 

1. Depth, 

2. Azimuth, 

3. Elevation, 

30 4. Radial component of blood velocity, 
5. Time. 
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Since a device of the present invention can automatically locate and lock onto the 
point in three dimensions having the maximum high-Doppler energy, i.e., maximum 
volume of blood having a significant radial velocity, unattended continuous blood 
velocity monitoring is one of its uses. By using the precise relative location of the 
5 point at which lock occurs as a function of depth, a device of the present invention 
can map the network of blood vessels as a 3-dimensional track without the hardware 
and computational complexity required to form a conventional ultrasound image. 
Using the radial component of velocity along with the three-dimensional blood path, a 
device of the present invention can directly compute parameters of blood flow, such 
10 as vector velocity, blood flow volume, and Doppler spectral distribution. 

A device having applications in a method of the present invention is a non- 
mechanical Doppler ultrasound-imaging sensor comprising probes, processing 
electronics, and display. Specific choices of probes allow the system to be used for 
transcranial Doppler (TCD), cardiac, dialysis, and other applications. 

15 Just as with other embodiments of the present invention set forth above, this 

embodiment has application for medical evaluation and monitoring multiple locations 
in the body. However, the transcranial Doppler application will be used as an 
nonlimiting example. Fig. 13 shows the overall TCD configuration and a typical 
definition of the display screen. The TCD system is comprised of one or two probes 

20 that may be attached to the head with a "telephone operator's band" or a Velcro 
strap. The interface and processing electronics is contained within a small sized 
computer. A thin cable containing from 52 to 120 micro coax cables, depending on 
the example probe design used, attaches the probe to the electronics in the 
computer. When the operator positions the probe on the head and activates the 

25 system, the Anterior, Middle and Posterior Cerebral Arteries and the Circle of Willis 
are mapped on the screen along with other blood vessels. The arteries or vessels of 
interest are selected by manually locating a cursor overlaid on the vessel 3-D map. 
The system locks onto the blood vessels and tracks their position electronically. A 
variety of selected parameters are displayed on the screen; e.g., the velocity, the 

30 pulse rate, depth of region imaged, gain and power level. Using only one probe the 
TCD can monitor multiple arteries (vessels) at a time. By way of example, presented 
on the screen are dual traces, one for each artery selected. The blood velocity can 
be dynamically monitored. As shown in Fig. 13 both the current blood velocity (dark 
traces) and any historic trace (lighter color) can be displayed simultaneously. The 

35 average blood velocity or estimated average flow for each artery is displayed below 
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the respective velocity trace. The image shows the arteries and the channel used for 
each artery. When two probes are used, the display is split showing signals from 
both of them. For example, using a different probe (i.e., different size) with the same 
electronics and display, the unit can be used to measure and monitor the blood flow 
5 in a carotid artery. Similarly, it can be used to perform this function for dialysis, 
anesthesia, and in other procedures. 

The sensor is a two dimensional array of transducer elements (e.g., piezoelectric 
crystals) that are electronically activated in both transmit and receive to effect a scan. 
For example, if a square (Wx/V) array is used, up to A/ 2 elements could be used at the 

10 same time. This is illustrated in FIG. 14 for the case of A/=8. The array need not be 
square. Any MxN array may be utilized in this manner. All received signals (52 in 
the example of FIG. 13) are sampled, digitized, and processed. This can be done, 
for example, in a desk top or lap top personal computer with additional cards for 
electronics and real-time signal processing as illustrated in FIG. 13 and FIG. 21. The 

15 array is phase steered or time-delay steered, depending on the bandwidth utilized, 
which depends in turn on the desired range resolution. The angular field of view 
shown in FIG. 15 is limited by the existence of grating lobes caused by array thinning 
(spacing the array elements more than wavelength apart). The concept is 
illustrated below for a 1 -dimensional array forming a beam that measures only one 

20 angle. For a two-dimensional array, this represents a horizontal or vertical cut 
through the cluster of beams shown in FIG. 15. 

The frequency utilized for TCD is usually at or near 2 MHz because higher 
frequencies do not propagate well through bone and lower frequencies do not 
provide adequate reflection from the blood cells. However, other frequencies have 

25 applications when examining other parts of the body. With a propagation velocity of 
1 .54 millimeters per microsecond, the wavelength is 0.77 mm. If a filled array is 
utilized, the element size and array pitch would be cf=0.77/2. For a cross-range 
resolution of 5.8 mm or less at a depth of 60 mm, the array size, L, must be at least 8 
mm (Resolution = depth x wavelength / L). Since N=Ud in FIG. 2, N must exceed 21 

30 and hence the array must have on the order of W 2 or over 400 elements. If the 
desired resolution is halved, the array size doubles and the number of elements 
exceeds 1,600. The array in FIG. 14 is said to be "thinned" because it only has 52 
elements. 

As explained above, "grating lobes" are ambiguities or extra, unwanted, beams 
35 caused by using a transducer array whose elements are too large and hence too far 
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apart. The following analysis illustrates grating lobe suppression for the worst case 
of narrowband signals and phase-shift beam processing. Time delay processing 
using wideband signals would be similar, but would further attenuate or eliminate 
grating lobes, resulting in even better performance. Naturally, one of ordinary skill in 
5 the art can readily program a processor to suppress or limit grating lobes with the 
equations described herein using routine programming techniques. 

Let 

x = (d/X)sme, (9) 
represent a normalization for the angle, ft, from which reflected acoustic energy 

10 arrives. The azimuth (or elevation) angle, ft, is zero in the broadside direction, 
perpendicular to the transducer array and d is the width (or length) of a single 
element of the receiver array. The wavelength of the radiated acoustic wave is X = 
c/f, where c is the acoustic propagation velocity (1540 meters/second in soft tissue) 
and /is the acoustic frequency (usually between 1 and 10 megahertz). The wide 

15 pattern in FIG. 16a is the element pattern 

a e (x) — sin tdc / 7ix . (10) 
The pattern is the product of the element pattern, the array pattern, and cos ft 

20 

a(0)=cos(0)a e {x)a a (x) (11) 

Each of the two component patterns is plotted separately as a function of ft in FIG. 
16a and the total pattern of equation (11) is plotted in FIG 16b. In the far-field, i.e., 
25 for Ar»L 2 t where r is the range or depth and L is the length of the aperture, the array 
pattern steered to the angle ft= ft is 

fl.W^^e^, (12) 

n«=0 

30 where w n is a weighting to reduce sidelobes and N is the number of elements in one 
dimension. As seen in Figure 16a, equation (12) is periodic in x. The peak at 9=6 0 
(6 0 = 0 in Fig 16) is the desired beam and the others are grating lobes. 

In the near field, when focused at (r^ft), equation (12) is replaced by the slightly 
35 better general Fresnel approximation: 
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>2/r «(w 0 )+ 




(13) 



10 



15 



20 



25 



30 



(provided that that the range significantly exceeds the array size, r>L), where x = d sin 
0/ A, as before, and 



Because the receiver aperture is sampled with a spatial period of d, the receiver 
array pattern will be periodic in sin 0, with a period of A/d (equation 12). This 
periodicity means that the array pattern is ambiguous. When the array is pointed 
broadside (0= 0), it will also be pointed at the angle 0= sin" 1 (A/d), for example. In 
terms of the normalized variable, x, the period is unity. Since | sin 0 \ cannot exceed 
1 , the variable x is confined to the interval [-d/A> d/A]. The conventional element 
spacing is d = All. Thus, in a conventional phased array, x is always between -0.5 
and +0.5, and hence ambiguities are not encountered. In a highly thinned array (d> 
X), there will normally be ambiguities or grating lobes as illustrated in FIG. 16a. The 
second grating lobe, at x=2 or 0= sin* 1 (2 A/d), is not real when d does not exceed 2A. 

FIG. 16b shows that the unsteered total pattern does not exhibit grating lobes. In a 
two dimensional array, the elements could be rectangular instead of square (d x x d y ), 
and the results would still be valid. Similar results could be obtained for an array in 
which the elements are staggered from row to row (and/or column to column). 

In FIG. 17 the same array is used as in FIG. 16, but the receiver element signals are 
combined with a phase taper that steers the beam to x = 0.2 or 6= 4.71°. In FIG. 
17b, we see that the grating lobes are not completely suppressed, with the largest 
one at x = -1+.2 = -.8 or 0= -19.18°. FIG. 18 shows this in decibels. The worst-case 
grating lobe is attenuated by at least 12 dB, even in the stressing case of extremely 
narrow band operation. These Figures were produced in MATLAB, using the 
following software (m-file): 

% M PATTERN mpattern.m Script to plot monostatic patterns vs. theta 

Mt=90; wave_length =0.77 ; d=1.875 , N=8, 

k=d/wave_length 

t = -Mt:0.1:Mt; 

tr = pi.*t./180; 

x=k*sin(tr) ; 

p=pi*x + eps; R=sin(p)./p; 



z = £ cos 2 



0/Ar. 



(14) 
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R=R. *cos (tr) ; 
n-0:N-l; 
% xo=0; 

xo=0.2; I steered 
5 e=exp( j*n' *2*pi* (x-xo) ) ; 
% w=hanning (N) ; 
% E=(2/N)*w'*e; 
E=(l/N)*ones(l / N)*e; 

subplot (211) ; plot(t, tabs(R);abs(E)]); 
10 ER=abs (E) . *abs (R) ; % Monostatic 

subplot (212) ; plot ( t, (abs(ER))); 
figure (2) ; plot (t, 20*logl0 (abs (ER) ) ) ; 
zoom on; 

15 The values of d and X used in the above example are representative for a 

transcranial Doppler application of the invention. If /= 2 MHz is chosen for the 
center frequency, the wavelength is 0.77 mm. An 8x8 array with a width and/or 
length of Z=15 mm, provides a one dimensional thinning ratio of 2 d/k = 4.87. A 15 
mm square array with half-wavelength elements would require more than 15,000 

20 elements. By thinning, this number was reduced to 52 provided that the angular field 
of view is limited to 2x4.71=9.42°. For a 1 cm array at 2 MHz, the hyperfocal 
distance (where the 3 dB focal region extends to infinity) is l}lAX = 3.25 cm. For a 15 
mm array, the hyperfocal distance is 7.3 cm. Thus, a fixed focus probe suffices for 
this application, but the quadratic phase distribution across the elements required to 

25 focus in depth should be added to the linear phase distributions required to steer the 
beams. 

Using the configuration described above, the cluster of beams in Fig. 15 is used to 
approximately locate the desired point for collecting the blood velocity signal. This is 

30 done initially, and is repeated periodically, in "mapping dwells" that are interspersed 
with normal measurement dwells. For example the output of each beam in the 
cluster would be Doppler processed by performing an FFT or equivalent 
transformation on a sequence of pulse returns. The pulse repetition frequency (PRF) 
would typically be less than or equal to 9 kHz to unambiguously achieve a depth of 

35 8.5 cm for the TCD application. In order to obtain a velocity resolution finer than Av = 
2 cm per second (to distinguish brain death), a dwell of duration as long as r= AJ(2 
Av) = 20 ms, or 170 pulses at 8.5 kHz, may be desired in the measurement mode. 
During monostatic mapping, 21 beams are scanned. If a mapping dwell is to be 
completed in 20 ms, only 8 pulses per beam are available, and an 8-pulse FFT would 

40 be utilized for each beam position. 
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The example shown in Figs. 16 through 18 was an 8 by 8 receiver array forming a 5 
by 5 cluster of beams. This is an example of an approximate rule of thumb for this 
invention, that an //element linear array is recommended for use in producing NI2 +1 
beams for N even and [W+1]/2 beams for Wodd. Thus, a 16 by 10 element 
5 rectangular array would preferably be used to form a 9 by 6 cluster of beams, though 
the actual number of beams formed is arbitrary. 

Because receive beams are formed only in a limited angular region, a wide-angle 
receiver element pattern (which usually implies a small element) is not required. In 

10 fact, the size of the receiver element can be as large as the element spacing. Thus 
the receiver array is "thinned" only in the sense that the element spacing exceeds a 
half wavelength. Since the element size also exceeds a half wavelength, the array 
area is not reduced. It is thinned only in terms of number of elements, not in terms of 
receiver area. Consequently, there is no reduction in signal-to-noise ratio, nor a 

15 requirement for increased transmitter power. 

FIG. 19 illustrates a means for increasing the angular field of view in the azimuth 
direction by extending the array horizontally. A similar scheme could be used 
vertically to extend the elevation F.O.V. The 52-element array of FIG. 14 becomes a 

20 single panel of the extended array. One panel is active at a time in FIG. 19. The 
beamwidth for FIG. 14, in radians, is nominally given by A/L At a range or depth of 
R, the cross range resolution is RM. (typically 3 to 5 mm). The F.O.V in millimeters 
at that same range is less than A//2+1 = 5 times that beamwidth. If a second panel is 
used in a planar configuration, the second panel translates the beam pattern to the 

25 right (or left) by the width of the panel, L=L 2 /2 (typically 8 mm). The field of view can 
be extended by more than this (it can even be doubled) by tilting the two panels in 
opposite directions to minimize the overlap in coverage of the two panels. 

Figure 19, with L, « L 2l simultaneously provides: (1) a large F.O.V. in the L 2 direction 
30 to allow for the simultaneous monitoring of two blood vessels more than an inch 
apart, (2) a large active array area for high sensitivity, and (3) a number of active 
elements below 60 and a total number of elements below 120. An alternative, shown 
in FIG. 20, has the array on the surface of a segment of a cylinder. This uses 52 
elements at a time with a total of only 84 elements (and hence only 84 cables). The 
35 Ln x L.2 active array translates around the curved surface as the beam is scanned 



WO 00/72756 



PCT/US00/14691 



65 

horizontally, if a symmetric F.O.V. extension (azimuth and elevation) is desired, a 
spherical surface could be utilized. 

Figure 21 is an overall block diagram depiction of the overall blood flow monitor. 
5 Most functions are performed by means of software in the digital processor. 

Naturally, one of ordinary skill in the art can readily program the processor to perform 
functions described herein using equations set forth herein and routine programming 
techniques. A possible implementation of the analog processing is diagrammed in 
FIG. 22. The A/D converter can be a bank of converters or one or more converters 

10 multiplexed amongst the 52 channels. If an extended array such as shown in FIG. 
19 or 20 were used, a switch would be included between the 52 processing channels 
in Figure 22 and the actual elements. Note that the 52 element array of FIG. 14 
represents an 8x8 array with corners removed (52=8x8-4x3). Other possibilities 
include a 24 element array (24 = 6 x6 -4x3). a 120 element array (120 = 12x12-4x6), 

15 etc. 

The transmitter produces pulses for each active element at a pulse repetition 
frequency (PRF) of 8,500 pulses per second. Each pulse will be at a frequency of f 0 
= 2 MHz and will have a bandwidth, B, of at least 250 kHz (e.g., a pulse no more 
than 4 microseconds long). 

20 

For measurement, only one or two beam positions need be insonated by a single 
probe. For mapping, many beam positions must be insonated, with several pulses 
on each for moving target indication (MTI) and/or Doppler processing. A 
measurement frame duration longer than 20 milliseconds (170 pulses at an 8.5 kHz 

25 PRF) may not be necessary because of the non-stationary (pulsed) nature of human 
blood flow. Mapping, requires several (4 to 11) pulses per beam position and many 
(e.g., 21 to 36) beam positions per frame. Since the Doppler resolution for mapping 
is not as fine as in the measurement mode, longer mapping frames can be used. If 
only 21 beams are formed with 8 pulses on each or if up to 34 beams are formed 

30 with only 5 pulses on each, a frame duration of 20 ms can be maintained even 
during search and mapping. 

FIG. 22 shows 52 identical receiver chains comprising 

1 . Processor controlled time gain control and time gate (open for up to 26 
35 microseconds for each pulse). 
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2. A limiter for dynamic range control. 

3. A low noise amplifier (LNA). 

4. A low pass filter (to pass \f\< f 0 +B/2 (e.g., |/) < 2.125 kHz) and reject \f\ > 5.875 
kHz by at least 40 dB (assuming f D =2MHz and 6=250 kHz). 

5 

AID conversion (typically 12 to 16 bits) is performed at an 8 MHz rate for each 
channel in FIG. 22. This keeps the analog filtering requirements extremely simple. It 
also permits extremely large bandwidths (up to 2 MHz) and time-delay steering. For 
narrower bandwidths and phase-shift steering, bandpass analog filtering and much 

10 lower sampling rates (determined by 8 rather than f 0 ) could be used. For the 8 MHz 
sampling rate, either time-delay or phase-shift beam steering can be utilized 
(depending on signal bandwidth). FIG. 22 depicts time delay steering for the 
transmitter. The distance from each array element to each focal point (each beam 
center at a nominal depth (e.g., 60 mm for TCD) would be pre-computed and stored 

15 either as a time delay or as a phase shift (depending on the type of steering) for each 
element for each beam. If phase shift steering were utilized on transmit, the 
transmitted signal could be created digitally in the processor, followed by D/A 
conversion for each element. Hence Figure 22 represents only one possible 
embodiment of the invention. 

20 

An example of the digital receiver processing for the case of an 8 MHz sampling rate 
per channel is described below. The input is 208 12 or 16 bit samples per pulse (8 
samples per microsecond x 26 microseconds to allow for a 4 cm deep radial mapping 
field of view), 8,500 pulses / second, and 52 channels. This results in a maximum 
25 average rate of 52x208x8500 = 91 .9 MegaSamples per second (or 1 .84 million 
samples in a 20 ms frame). During measurement, the range interval can be 
narrowed to less than 1 cm, reducing the number of samples per pulse to 32. The 
average rate for measurement and monitoring becomes 14 megasamples per 
second. The receiver processing steps are as follows: 

30 • Buffer (to allow subsequent processing to be performed at the average rate). 
• Digitally Down Convert to Baseband (make I and Q). 52 channels in parallel. 
Multiply input samples by samples of a 2 MHz cosine wave and -sine wave to 
create In-phase and Quadrarture samples, respectively. Since the samples are 
% cycle apart, the multiplicands are all 0, 1, or -1, and hence no multiplications 
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are needed. If r(j,p) is the real p m sample from the f channel, the complex low-pass 
signal, s(j, p), has a real part for p = 0,1, 2, 3, 4, 5,... given by 

rftO), 0, -/t/,2), 0,/V i 4) i 0 i ... 
and an imaginary part given by 

0, -rO',1), 0.fft3), 0, -ftf.5).... 
This provides a data rate 2 times the input rate because the data is now complex. 

• Pre-Decimation Low-Pass Digital Filter. Filter 52 complex channels. 

Pass |/| < S/2, reject \f\ > r- B/2, where r is the sampling rate after sample rate 
decimation (e.g., 1 MHz). If B = 250 kHz, r could be as low as 500 kHz. If B is 
large, r could be 2 or 3 MHz. If receiver phase-shift steering were to be 
performed, the output samples would be computed at the decimated rate. If 
receiver time-delay steering is to be used, we output 8 million complex samples 
per second and postpone sample rate decimation until after beam formation. 

• Perform MTI or create coarse Doppier cells. For every channel and every 
range sample, either digitally high-pass filter the sequence of pulse returns to 
suppress clutter from tissue and bone or perform 52x208 8-point discrete Fourier 
transforms (DFTs or FFTs) for each mapping frame. (Six points of the 8-point 
complex DFT provides 3 positive and 3 negative coarse Doppier cells.) 

• Perform Digital Beamforming. Case 1: Time Delay Beamforming with Sample 
Rate Decimation uses a set of pre-computed time delays to reduce 52 complex 
channels with 208 samples per pulse to one of M (e.g. 21) complex beam outputs 
with 25 samples (range cells) per pulse. The example given here assumes 8:1 
decimation. 

The maximum delay is slightly less than 0.75 us = 6 r, where T= 1/8 microsecond 
is the time between input samples. For a given pulse return, the sample (*= 1, 
2, 25) of the i* beam, i = 1, 2, .... M, is denoted by b(i, k). The /? m sample (p = 
1, 2, ... f 208) of the j th input channel (/ = 1, 2, 52) is denoted by s(j,p). Let<% 
be the delay required for the signal in channel j to produce beam i. 
For a given pulse return, the k* complex 1 MHz rate output sample for beam i is 

b(i,k) = f>, s{j\ &[k + 1] - b g ) + (1 - a 9 )sO\ ** + !]- b y -1)} 
where by is the integer part of dy/T (between 0 and 6) and a i} is the fractional 
remainder (between 0 and 1). Determine power or amplitude in each output 
Doppier bin as I 2 + Q 2 or its square root : 
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Case 2: Phase-shift beamforming of already decimated data involves only a 
sequence of inner products of 52-dimensional complex vectors of element 
values with a complex vector of representing the required phase shifts. 

• Display Coarse Blood-Vessel Color- Flow Map. Coarse blood vessel map is 
5 the set of range, azimuth, and elevation cells with high power, with 6 Doppler 

values. Blue and red represent positive and negative Doppler, with saturation 
related to radial velocity and intensity related to power. 

• Initialize Acquisition. The user, looking at an azimuth-elevation Coarse Map 
(with depth automatically truncated to a set of values that should include the 

10 MCA), moves the transducer and looks for a high-intensity, saturated spot. He 
can center the probe on that spot or he can have a device of the present 
invention display a range interval corresponding to the ACA, in which case he 
can make sure that both vessels are well within the angular field of view of the 
probe. 

15 • Acquisition and Tracking of one or two points being monitored. This is done 
with a single transmit beam focused on the spot identified above for several 
frames. Digital Down-conversion, low-pass filtering, and MTI are performed as 
before, but beamforming is different. Five receive beams are simultaneously 
formed. These are a sum beam and four monopulse difference beams, all 

20 steered to the same point as the transmit beam. Each monopulse beam is 

equivalent to the difference between the outputs of a pair of beams displaced on 
opposite sides of the focal point. The four monopulse pairs are in 45 degree 
intervals with the first being horizontal, and the third being vertical. The 
monopulse-difference output with the largest magnitude is divided by the output 

25 of the sum beam. The imaginary part is the "monopulse ratio" used to re-steer 
the beam (in the difference pair direction) so that it is better centered on the 
vessel. This procedure can be repeated in an effort to drive all four monopulse 
ratios to zero. 

• Measurement and Tracking. Tracking continues as described above during the 
30 measurement mode. Measurement is made with fine Doppler resolution (128 

point FFT) applied to only the sum beam. In a 15 ms frame, data from 128 
pulses are collected (52 channels, 6 range samples). The pulses are Hamming 
weighted and FFTd. This produces 128 Doppler bins (for each range bin and 
element), 66.4 times a second. Real sum beam outputs would then be produced 
35 (using monopulse-guided steering) for each of 64 to 126 of these Doppler bins. 
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• Track maintenance and re-acquisition. Tracking is continued in parallel with 
measurement. If a monopulse ratio suddenly deviates far from zero and is not 
brought back to zero in one or two iterations, loss of track is declared. Re- 
acquisition is attempted autonomously by re-steering the beam by an amount 
5 determined by correlating a current color flow map with a stored earlier version 
(from before track was lost). If this is unsuccessful, (monopulse ratios do not all 
converge to zero) an alarm is sounded so that the user can return to repeat 
initialization of acquisition. 
Correlation with previous maps will be periodically applied to prevent wandering of 
10 the data collection point along the vessel being tracked. 

For tracking purposes, a monopulse tracking method described above can be used. 

FIG. 23 illustrates the segment of a vessel in a single resolution cell, after rotation. 

15 The resolution cell shown is not a cube because the range resolution might be finer 
than the cross-range resolution. The illustrated circular cylinder represents blood 
cells in a vessel reflecting energy at a fixed Doppler frequency. These represent a 
cylindrical annulus of blood cells, at a constant distance from the vessel wall, moving 
with approximately the same velocity. In the single resolution cell of Fig. 23, the 

20 return at the highest Doppler would represent a line in three-dimensional space (the 
axis of the vessel) and hence a point on the azimuth axis after rotation. When 
applied to the highest Doppler output, the Sum beam would have broad peak at zero 
azimuth (a * 0) and the monopulse ratio, r=Az/Sum, will be a linear function of the 
azimuth angle to which the array is phase steered: 

25 r(a)=ka. 

This result can be attained by applying the same phase across the aperture for the Az 
and Sum beams, but using the derivative of the Sum beam amplitude weights with 
respect to x and y respectively for the Az and El aperture weights. 

30 Many other variations and modifications of the invention will be apparent to those 
skilled in the art without departing from the spirit and scope of the invention. The 
above-described embodiments are, therefore, intended to be merely exemplary, and 
all such variations and modifications are included within the scope of the invention as 
defined in the appended claims. 



